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We study the effect of quantum fluctuations on the half-polarized magnetization plateau of a 
pyrochlore antiferromagnet. We argue that an expansion around the easy axis limit is appropriate 
for discussing the ground state selection amongst the classically degenerate manifold of collinear 
states with a 3:1 ratio of spins parallel/anti-parallel to the magnetization axis. A general approach 
to the necessary degenerate perturbation theory is presented, and an effective quantum dimer model 
within this degenerate manifold is derived for arbitrary spin s. We also generalize the existing 
semiclassical analysis of Hizi and Henley [Phys. Rev. B 73, 054403 (2006)] to the easy axis limit, 
and show that both approaches agree at large s. We show that under rather general conditions, 
the first non-constant terms in the effective Hamiltonian for s > 1 occur only at sixth order in 
the transverse exchange coupling. For s > 3/2, the effective Hamiltonian predicts a magnetically 
ordered state. For s < 1 more exotic possibilities may be realized, though an analytical solution of 
the resulting quantum dimer model is not possible. 

PACS numbers: 75.10.-b,75.10.Jm,75.25.+z 



I. INTRODUCTION 

Magnetism is an inherently quantum mechanical ef- 
fect that ultimately arises from the exchange processes 
in an interacting many-particle quantum system. Thus, 
magnetic systems can reveal much about the richness of 
quantum mechanics itself, especially when cooperative 
effects are operational and emergent effective low-energy 
theories describe the relevant physics. A few notable the- 
oretical examples are the spin liquid phases obtained in 
frustrated magnetsiaSi&i^. Quantum effects may play a 
central role in other types of order as well. Recent exper- 
iments on a number of insulating chromate compounds, 
namely CdCr204 and HgCr204, have shown peculiar 
features in the low-temperature magnetization as a func- 
tion of applied magnetic field. At low temperatures the 
magnetization grows linearly with magnetic field up to 
some critical value at which point there is a sharp jump 
in magnetization onto a rather wide plateau with half 
the full saturation magnetizationiSiLS With sufficiently 
large fields it is possible to observe a smooth transition 
off the half-magnetization plateau and a gradual increase 
in magnetization up to what may be a fully polarized 
plateau stated As described in Ref. 0, it is expected that 
the magnetism in these compounds is well described by 
the Heisenberg antiferromagnet (AFM) of spin s = | on 
the pyrochlore lattice. 

Motivated by these experimental examples of a half 
polarization plateau in the pyrochlore Heisenberg AFM, 
we conduct a theoretical study of the quantum py- 
rochlore Heisenberg AFM for any spin value s, in a strong 
magnetic field, focusing on a half polarization plateau. 
While the physics of HgCr204 is probably determined 
to a large degree by the classical physics of spin-lattice 
interactions;ifi^i*i2iiii4 it may be that in other similar 
compounds where coupling to phonons is weak, quantum 



effects could play a significant role. 

In any case, the general problem of determining 
the spin state on plateaus of non-zero magnetization 
in frustrated magnets occurs in a large number of 
materialo 15 i 16 i 17 i 18 i 19 i 20 i 21 i 22 . At fields large enough to 
induce substantial magnetization, the ground state is 
expected to be very different from the zero field state 
and one would ideally pursue a theoretical approach that 
takes advantage of the large external field. The methods 
developed indeed use this explicitly, and the particular 
application to the half-polarized pyrochlore magnetiza- 
tion plateau provides a rather non-trivial test bed. We 
make use of the large field to justify an easy-axis approx- 
imation to a nearest-neighbor XXZ antiferromagnetic in 
an external field. Physically, at large fields the spin is 
oriented on average more along the field axis than trans- 
verse to it. Furthermore, specifically on a magnetization 
plateau, general arguments imply that the static trans- 



verse moment vanishes on every site, (S^ 



0. Thus 



we expect that Degenerate Perturbation Theory (DPT) 
about the easy-axis (Ising) limit should be justified on 
the plateau, and one may thereby derive an effective 
Hamiltonian. This effective Hamiltonian acts in the con- 
strained "3:1" space of states with 3 majority spins with 
Sf = +s and 1 minority spin with Sf — —s on each tetra- 
hedron. This space is macroscopically degenerate, and all 
its members have half the saturation magnetization. 

The reader may well wonder whether there is any need 
for an approach of this sort, given the successes of the 
large s semiclassical spin-wave method in many other 
contexts. Indeed, for unfrustrated antiferromagnets, it is 
known that the 1/s expansion gives reasonably conver- 
gent results even down to s = 1/2. However, this conver- 
gence is strongly dependent upon the lattice - large cor- 
rections to the spin-wave dispersion have recently been 
obtained even for the rather weakly-frustrated triangu- 
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lar latticei^iSa In highly frustrated magnets such as the 
pyrochlore, another approach is warranted. Particu- 
larly worrisome in the 1/s expansion is the difficulty 
of treating tunneling, which is non-perturbative in this 
method2£i2§*2!. By contrast, in the easy-axis expansion, 
tunneling and virtual exchange are treated on the same 
footing. Of course, for large s both approaches must 
agree, and we will indeed check that this is the case in 
our specific application. 

The effective Hamiltonian describing splittings within 
the degenerate manifold of Ising ground states will gener- 
ally take the form of a constrained quantum Ising model. 
As explained in a previous publication^ for the case of 
the pyrochlore half-magnetization plateau in the easy- 
axis limit, it can be cast in the form of a "Quantum Dimer 
Model" (QDM) on the bipartite diamond lattice. Such 
QDMs are known to display both ordered and disordered 
(spin-liquid) ground states in different regions of their 
phase diagrams iSi^S2i22i2Si^i We derive the parameters of 
this QDM for general s, and discuss several limits and the 
expectations for the plateau ground state. For simplicity 
of presentation we perform this calculation here for the 
simplest XXZ spin model with no additional anisotropy 
or other spin interaction terms. However, the method 
is straightforwardly generalized to include other on-site 
(e.g. uniaxial anisotropy^Si) or nearest-neighbor (e.g. bi- 
quadratic) interactions without substantial increase in 
computational complexity. More generally, the flexibility 
to include such effects allows one to consider the quantum 
effects upon the ground state selection within a magne- 
tization plateau even when the dominant mechanism of 
plateau stabilization is a classical one. 

A remarkable feature of the DPT is that all diago- 
nal terms describing splitting of the low energy manifold 
vanish below sixth orderl For s > 1, off-diagonal tunnel- 
ing terms also vanish up to this order, so that the entire 
effective Hamiltonian is determined by terms of sixth or- 
der and higher. This behavior is similar to a result of 
HenleySS, that in the large- s limit, the effective Hamilto- 
nian is expressed entirely in terms of a "spin flux" vari- 
able involving a product of 6 spins around an elementary 
loop of the lattice. We show here that our result is rather 
general, and originates from two basic features: the ab- 
sence of non-trivial loops of length less than 6 links on the 
pyrochlore lattice, and the fact that all low-energy spin 
states on a single tetrahedron are permutations of one 
another. From our proof of this result, it can readily be 
seen that similar behavior holds for any lattice of corner 
sharing simplexes with only on-site and nearest-neighbor 
interactions and permutation-related ground states on a 
single simplex. We will apply the methods of this paper 
to other such problems of interest in future work. 

For the pyrochlore magnetization plateau and QDM 
studied here, the conclusions are as follows. For s > 3/2, 
we find that diagonal terms in the QDM are much larger 
than off-diagonal ones. In this case, the latter are negli- 
gible, and because the diagonal QDM is effectively clas- 
sical, it is soluble and the ground state is ordered. We 



discuss the preferred spin ordered states as a function of 
s. For s < 3/2, the off-diagonal terms are non-negligible, 
and a simple solution is no longer availed. For s = 3/2, 
various arguments lead us to still expect an ordered state. 
After correction of an error in Ref|9|, we find two candi- 
date states for this case. One of these is the R state 
discussed previously in Ref0, which is the state contain- 
ing the maximal number of hexagonal loops with alter- 
nating spins (flippable plaquettes in the QDM language). 
Another candidate is a -\/3x V3 state with a planar struc- 
ture. In fact, the diagonal terms in the effective Hamil- 
tonian do not entirely fix the relation between adjacent 
planes in the latter state, so there is additional degen- 
eracy whose breaking we cannot resolve at the present 
time. Because the off-diagonal and diagonal terms are 
comparable in this case, however, some other states may 
also be possible, and a definite conclusion must await 
more serious computational (e.g. quantum Monte Carlo) 
analysis. For s < 1, the off-diagonal term in the QDM is 
dominant. In this case, either the R state or U(l) spin 
liquid is the most likely candidate ground state. Indeed, 
as argued in Rcf 33, it is quite possible that the simplest 
QDM displays a direct quantum phase transition between 
these two states. 

The ground state of the QDM just discussed is de- 
termined only by the dimensionless ratios of coupling 
constants. However, the DPT calculation also gives the 
overall scale of the effective interaction in terms of the 
microscopic exchange J. For s = 3/2, the largest in- 
teraction energy (extrapolated from the easy axis per- 
turbation theory to the Heisenberg limit) generated by 
quantum fluctuations is only w 0.02 J. Were this the true 
scale for ground state selection in the degenerate 3:1 man- 
ifold in HgCr204, the magnetic ordering would occur at a 
temperature of this order, i.e. ss 0.2K. Experimentally, 
however, magnetic ordering is observed at a substantial 
fraction of the temperature of onset of the plateau forma- 
tion, which is around 6K. The closeness of the ordering 
and plateau scales in experiment suggests that both are 
determined by the same physical mechanism, and argues 
against the importance of quantum fluctuations in the 
ground state selection in HgCr 2 04. Indeed, we have re- 
cently shown 10 that the same spin-lattice coupling which 
leads to plateau formation can also account for the state 
selection. Curiously, the R state is also stabilized by 
the lattice mechanism. This is symptomatic of the very 
strong constraints defining the 3:1 QDM states, which 
lead rather different microscopic interactions to favor the 
same ground state. For s — 1 and s = 1/2, the DPT 
gives much larger characteristic scales for the QDM, the 
off-diagonal term being of order 0.16J and 1.5J in the 
two cases. Thus such s < 1 antiferromagnets, if realized 
experimentally, would be promising systems to observe 
quantum fluctuation effects. 

This manuscript is organized as follows. In Section ITU 
we describe our theoretical model, the nearest-neighbor 
quantum Heisenberg antiferromagnet on a pyrochlore lat- 
tice in an external field. An easy-axis limit is taken un- 
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der the assumption of the suppression of transverse spin 
fluctuations in large magnetic fields. After applying de- 
generate perturbation theory (DPT) in the transverse 
spin fluctuations, an effective dimer model emerges in 
Section Ifffl that can be used to obtain an approximate 
ground state of the original model. In Section wc 
carry out a large s analysis of the XXZ model, deriving 
a different effective Hamiltonian splitting the 3:1 mani- 
fold of degenerate states. This new effective Hamiltonian 
turns out to coincide with the s — ► oo limit of the effec- 
tive Hamiltonian from the DPT analysis. In SectionlVlwe 
explore the ground state of the diagonal part of the effec- 
tive Hamiltonian from DPT. In Section fV Bl we explore in 
more generality the appropriate Quantum Dimer Model 
(QDM) of which all our effective Hamiltonians are spe- 
cial cases. We conclude the main text of this paper with 
a discussion of our results in Section IVII In appendix 1X1 
we analyze how the half polarization plateau is modi- 
fied by quantum fluctuations. An alternative method of 
performing DPT is presented in Appendix 151 and shows 
perfect agreement with the result of Section UTTl Finally, 
in Appendix IO we explore the states degenerate with the 
y/3 x y/3 states, found for s = |. 

II. MODELS 

A. Hamiltonians and Limits 

We begin with the simple spin-s Heisenberg antifer- 
romagnet (AFM) residing on the sites of the pyrochlore 
lattice in the presence of a magnetic field H, 

(ij) 3 

On the pyrochlore lattice one may recast the nearest- 
neighbor exchange in terms of the total spin on tetrahe- 
dra using the identity 

2]>>.S,=5» 2 -5:E( S J ) 2 ' (2) 

(ij) t t jet 

where St = Ylj^t ^j 1S ^ ne sum °f spins on a tetrahedron 
labeled by t, and {Sj) 2 \j) = S(S + l)\j). This gives the 
more convenient form 

n = { E h ) 2 - h2 ] - ( 3 ) 

( 

where we have introduced the dimensionless magnetic 
field h = H/2J = hi, and ignored a trivial constant 
term in the Hamiltonian. 



1. Classical limit 

The form in Eq.© makes the behavior in the large s 
limit apparent. In this limit the spins behave classically, 



and one may replace Si — > shi, where hi is a unit vector. 
The ground states then consist simply of all states for 
which St = s^2 iet hi — h on every tetrahedron. This set 
has a large continuous degeneracy. Furthermore, since 
the magnetization is simply half the sum of the St (be- 
cause each spin is contained in two tetrahedra), this im- 
plies a continuous linear behavior of the magnetization 
with field. Thus, in this model magnetization plateaus 
can emerge only from quantum corrections to the classi- 
cal limit. 



2. Easy axis limit 

An alternative approach exploits the fact that, with 
the application of the magnetic field, the global SU (2) 
symmetry of the bare Heisenberg model is broken down 
to a U(l) symmetry (rotations about the magnetic field 
direction). Moreover, when the magnetization per spin 
is substantial, on average the transverse components & 
are smaller in magnitude than the longitudinal ones. It 
is therefore natural to treat transverse and longitudinal 
exchange couplings on a different footing, with the latter 
taking the dominant role. Formally, this is accomplished 
by replacing the isotropic Heisenberg Hamiltonian by an 
XXZ model: 

H = Ho+H ll (4) 

where 

Wo = yE [( s ? ~ h ) 2 ~ h *\ ~ J *E (^) 2 • ( 5 ) 
t i 

and 

«i4e( s ^+^) ■ ( 6 ) 

(ij) 

We use the notation S% — J2i<£t &i > an d we nave made 
use of the identity 

E^ = ^E(^) 2 -E( 5 ?) 2 ■ ( ? ) 

(ij) t i 

In the equations above, and elsewhere in this manuscript, 
(ij) denotes a sum over nearest neighbor sites on the py- 
rochlore lattice, and Sf are the spin ladder operators. 
Note that in the Heisenberg model J± = J z — J, but the 
more general XXZ model has all the same symmetries as 
the former even when this condition is not obeyed. From 
the above reasoning, we expect that the transverse terms 
involving Jj_ may be treated as "small" perturbations in 
the strong- field regime of interest. We note that this is 
expected to be a particularly good approximation when 
the system exhibits a magnetization plateau. This is be- 
cause, as described in the introduction, (S^) must vanish 
in such a state. Formally, this "easy axis" limit consists 
of taking J± <C J z and doing degenerate perturbation 
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theory in a = We will assess the validity of this 

approximation later by considering the magnitude the 
perturbative corrections extrapolated to a = 1. Finally, 
we note that several other effects can stabilize a collinear 
state. One is the addition of easy axis anisotropy, 



(8) 



with K 3> Jj_ > 0. A second mechanism is biquadratic 
exchange, 



(ij) 



(9) 



with b > 0. A term of this form can be generated dynami- 
cally from spin-lattice interactions, known to be strong in 
HgCr204. The DPT treatment discussed below can read- 
ily be generalized to include either or both of the terms in 
Eas. l|8l9|) . For simplicity of presentation we do not do so 
here. While easy-axis terms similar to Eq. (JHJl are allowed 
for s > 1/2, this particular simple form, with the same 
spatial direction for the local easy axis of all spins, is not 
physically appropriate for the cubic pyrochlore spinels, 
and the proper anisotropy terms allowed by symmetry in 
these materials are likely to be very small in any case. 



B. Magnetization process in the Ising model 

The evolution of the ground state with field in the ex- 
treme easy-axis limit a = is less trivial than in the 
classical limit. The system is then described by the Ising 
Hamiltonian, Eq.JSJ). We shall focus on the h > case, 
as the case ft < is equivalent. The expression Eq.JSJ) 
can be written as a sum over tetrahedra Ho = Y2t TLt 
with 



OS? ~h) 2 -h 2 -J2 ( S j 



jet 



(10) 



and therefore if one can minimize this energy on each 
single tetrahedron, one will have attained the minimum 
energy of the many-body system. 

The magnetization S* of any individual tetrahedron 
is quantized to the values Sf — 0, ±1, ±2, ... ± 4s. The 
1st term in (llOfl favors the magnetization of the tetra- 
hedron to take on a value close to the integer part of 
the dimensionless magnetic field. The 2nd term in i|10|) 
favors larger z-components of the spin values S? = ±s. 
This 2nd term is trivial only in the spin i case where the 
spin-1/2 Pauli matrices square to the identity. 



A state with the magnetization S t z 



[ft] (the 



integer part of h) clearly minimizes the energy of the 
first term in (|10l) . However, given the magnetization, 
there is some freedom for the values of the spins on each 
tetrahedron. The 2nd term in (|10|) reduces this freedom 
by adding an energy cost for small S z components. 



Consider 4 spins on a tetrahedron with total mag- 
netization to, and individual values of Sf, Sf, Sf, Sf. 
Now compare the energy of this state with that of 
Sf — l,Sf + 1, S3, Sf . The only energy difference comes 
from the second term in (|10|) 



AE 



(si + iy 

2(1- Sf 



-(Si 
s 2 2 ) 



iy + (siy + (szy 



(11) 



From this expression one deduces that if one begins with 
Sf > Sf, then it is energetically favorable to increase 
Sf even more at the expense of Sf. This increase in 
Sf can only be halted by one of Sf 2 reaching an ex- 
treme spin value of ±s. From this reasoning we conclude 
that the lowest energy state on a single tetrahedron with 
a fixed total magnetization to has a spin configuration 
with the largest possible number of extreme valued spins. 
This also makes intuitive sense, since ideally all the spins 
should take on extreme values if possible. In three par- 
ticular choices of to, all the spins take on extreme values. 
Zero magnetization to = with s, s, —s, —s, half polar- 
ization of to = 2s with s, — s, — s, — s and full polarization 
to = 4s with s, s, s, s. For to < 2s we find the lowest 
energy configuration for the spins s, s, —Sjin — s, and for 
to > 2s, we find s, s, s, to — 3s. Finally, we can find the 
minimal energy for given magnetization to by using the 
spin configurations described above for every value of to, 
and plugging them into (|10|) . For to < 2s we find 



E m = ^ [2to(s - ft) - 4s 2 ] 



and for to > 2s we find 



E,, 



2 



[2(m-2s)(3s-ft)-4fts] 



(12) 



(13) 



From the above expressions it is easy to see that at ft = s 
all to < 2s yield the same energy - all these states are 
degenerate at this field value. Similarly, for ft = 3s all 
to > 2s yield the same energy. For other values of the 
magnetic field, we find for h < s the lowest energy state 
is the m = zero magnetization state, for ft < s < 3s 
the lowest energy state is the to = 2s half polarized state, 
and for ft > 3s the lowest energy state is the to = 4s fully 
polarized state. 

The three lowest energy states in the various magnetic 
field ranges have all spins at extreme values ±s, and can 
be realized on every tetrahedron in the pyrochlore lat- 
tice. The to = state induces a degenerate manifold of 
states with every tetrahedron having s, s,—s, —s on it. 
This 2:2 proportionality is well known as the "ice rules" 
encountered in a particular phase of water ice^ as well 
as spin ice compounds^ The half polarization states are 
also massively degenerate, with every tetrahedron in a 
3:1 proportionality of S| = +s to SJ = — s spins (or 3 up 
one down). This particular degenerate manifold will be 
the focus of the remainder of our discussion. To summa- 
rize, the magnetization curve for JSJ exhibits 3 plateaus 
at zero, half and full polarization, for all values of s. 
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III. EASY AXIS DEGENERATE 
PERTURBATION THEORY 

A. Structure of perturbation theory 

1. Basic formulation 

In the previous section, we observed that the extreme 
easy axis limit of the Heisenberg model exhibits a broad 
magnetization plateau at half polarization. However, the 
ground states on this plateau are macroscopically degen- 
erate, consisting of all states with a 3:1 ratio of majority 
and minority spins on each tetrahedron. In this section, 
we study the splitting of this degeneracy by perturba- 
tion theory in J±. We employ the following formulation 
of Degenerate Perturbation Theory (DPT). Define the 
projection operator, V, onto any degenerate manifold of 
states M. Consider any exact eigenstate |\&). Its pro- 
jection |$o) = V\$l) satisfies the "effective Schrodinger 
equation" 



Eq 



VHi Q nv 

n=0 



|*o) =£1*0) = Wcff|*o), (14) 



where the operator Q = E _ Ha (1 — V)Hi. Because the 
resolvent contains the exact energy E, Eq. l(T4"|l is actually 
a non-linear eigenvalue problem. However, to any given 
order of DPT, E may be expanded in a series in Jj_ to 
obtain an equation with a true Hamiltonian form within 
the degenerate manifold. Each factor of Q is at least of 
0(J±) due to the explicit factor in Hi, with higher order 
corrections coming from the expansion of E. Once |^o) 
and E are known, the full wavefunction can be recon- 
structed as = (1 - g)- 1 |* ) = EZo^o)- 

Considering the lowest order term in DPT that breaks 
the degeneracy, the precise energy E = Eq + 0(a) in the 
resolvent can be replaced by Eq, where 0(a) represents 
possible energy shifts from lower order terms that do not 
break the degeneracy, and Eq is the 0-th order energy of 
the degenerate manifold of states. 



2. Order of off-diagonal terms 

Every order in DPT can in principle have diagonal (in 
the Sf basis) as well as off-diagonal terms in which the 
degeneracy is removed. Any off-diagonal term in the ef- 
fective Hamiltonian must flip spins in such a way as to 
preserve the 3:1 constraint on each tetrahedron. This 
can only be accomplished by flipping spins around a 
non-trivial closed loop on the pyrochlore lattice (see, e.g. 
Ref 2). The smallest such loop involves flipping spins on 
3 different bonds, and flipping a spin from S z — ±s to 
S z = =Fs requires Hi to act 2s times, so off-diagonal pro- 
cesses occur first at order O(J^). Therefore, below this 
order of DPT, one need consider only diagonal terms. 
In subsection IIII Bl we will demonstrate that the lowest 



order diagonal energy splitting term for any s can oc- 
cur only at 6th order. For spin s — h an off-diagonal 
term appears at 3rd order in DPT, and no diagonal en- 
ergy splitting occurs at this order, resulting in a purely 
off-diagonal effective Hamiltonian. For spin s — 1 the 
lowest order diagonal and off-diagonal terms simultane- 
ously appear at 6th order. For any higher value of s, 
the diagonal energy splitting appears at a lower order 
than any off-diagonal term can occur, and therefore the 
leading order effective Hamiltonian is purely diagonal in 
the 3:1 states. We will nevertheless compute the first 
non- vanishing off-diagonal term for various values of s in 
Section IV Bl to use its magnitude for an assessment of 
the validity of the truncation of the DPT expansion. 



3. Unitarily transformed formalism for diagonal terms 

We next develop a scheme to compute the diago- 
nal terms, by unitarily transforming the expression in 
Ea.i|14|l to obtain a formula for the diagonal effective 
Hamiltonian with all dependence upon the spin state ex- 
plicit. The 3:1 manifold can be described using Ising 
variables to indicate which spins are minority sites. That 



djS with o~j 



±1 



is, in the 3:1 states, we denote S? 
the Ising variable. At n th order, presuming that all lower 
order terms are constants, the diagonal terms in the effec- 
tive Hamiltonian constitute the function of the set {<Ji} 
given by 

tun*}] = (-i) n+1 MW]l (Wifter^wiMfo}]), 

(15) 

where the resolvent 1Z = (Ho — Eq) 1 , Q = l — V, and 



M{«7i}]> = ®i\S? = <T iS ). 



(16) 



The assumption that all lower order terms are constant 
allows us to replace E by Eq in the denominators in 
Eq. (|14|) , since the constant corrections to E lead to higher 
order terms in the effective Hamiltonian. 

The dependence upon the <Xj in Eq. l|15|) . is not explicit, 
but, following Hizi et al^&, it can be made so by a unitary 
transformation. The operator 



(17) 



effects a rotation about the a;-axis in spin space only for 
the minority spins. This interchanges raising and lower- 
ing operators, and reverses the orientation of Sf for these 
sites. We may therefore write 



mWi}}) = urn 



where 



(18) 



(19) 



is the fully polarized state, which is now independent of 
o~i . Then we have 

= (-l)" +1 (^o| [HiHQ) Wihfo), (20) 
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where 



6 = U*OU, 



(21) 



for any operator O. In what follows, all the operators 
appearing in Eci. (|20|) above, will be simplified so that 
their dependence upon <Ji becomes explicit. 

First consider H%. It consists, from Eq.JBJ), of a sum of 
operators transferring spin 1 between two nearest neigh- 
bor sites, i.e. a bond of the pyrochlore lattice. We define 
the nearest-neighbor connectivity matrix of the lattice 
Tij — Tji = 1 when i and j are nearest neighbors, and 
Tij = otherwise. With this terminology we write © as 



(22) 



After the unitary transformation, one obtains 

(1 + <Xj<7 



(S+S- + h.c.) 



(23) 



Here the expressions ^ 1±a ^ a ^ denote "Ising delta func- 
tions" that select the cases in which the two a^j have 
the same or opposite sign. 

Assuming the lowest order term in DPT that splits the 
3:1 configurations is a diagonal term of order no, the only 
3:1 configuration which can be reached as an intermedi- 
ate state in Ea. (|15f) for any n < tiq is the starting state 
lV'[{ (7 i}])- Under the unitary transformation, this state 
maps to \ipo), and therefore the projection operator Q 
may be replaced by 



G->i-hfo)<ito| 



(24) 



in Ea. lj2U|) . 

Finally, we consider the resolvent. Using U^S*U — 
crjSf, one finds 

^ = y^ TijWjSiS'j - 2J z h VjSj - Eo- (25) 

First we note that because both Hq and Hi conserve the 
total magnetization of the lattice (this is just the con- 
served quantity arising from the global U(l) symmetry), 
the term J2j a j &j remains unchanged at every stage in a 
DPT process, and we can therefore absorb this term into 
the constant energy Eq. Clearly, the inverse resolvent 
should vanish when acting upon the fully polarized state 
|V>o)- Hence we may absorb the constant energy Eq into 
the sum as 



2 



(26) 



We can simplify the resolvent in the restricted space 
of virtual states which will be accessed in evaluating 
Eq.(j20|). In particular, the ct; configurations are re- 
stricted to the 3:1 manifold. Furthermore we note that 
all intermediate states will have only some small finite set 
of spins whose Sf quantum numbers arc different from 
s, due to the action of H\. Let us consider then the ac- 
tion of the resolvent on a state for which this set of sites 
is denoted by F. In this case, only terms in Ea. (|26f) for 
which at least one of i or j is in F can contribute. Thus 

^ = fE v m (sis* - s 2 ) 

One may replace the sum over j by X^f = 2~2j ~ SjeF 
to obtain 

K- 1 = y ]T T^a, (Sf - S ) (SJ - S ) 

ij'SF 

+ J z s g a t j (5? - s) . (28) 

The crucial observation is that the 3:1 constraint implies 



£rVj =4-2^ 



(29) 



This is because once Oi is specified, the set of its neigh- 
bors is also specified (see also Fig. El - Eq ■ (|29|) allows one 
to eliminate the latter sum and obtain 



n- 1 = ^^r iiC r^(5?- s )(5|- s ) 

ijGF 

+2J z s ^(2ai-l)(S? -s). 

JGF 



(30) 



Using again the observation that y\ a jSj remains un- 
changed throughout the stages of any DPT process, it 
is equal to the constant J2j a jSj = J2j a i s - Using this 
fact, we finally obtain 



£ -1 = ^E r ^^(^- s )(^- s ) 



ijSF 

-2J z sJ2(S? ~ «) 



(31) 



By successive action of Hi, Q, and 1Z using 
Eqs. (|23I24I31|1 . one can obtain explicit expressions for 
any intermediate state in the DPT expression of Eq. (|20|l , 
with n < Uq. For example, one action of each of these 
operators gives 



(32) 



as 



4(4j 1) ^~"^ aia2 ^ G a.\ a ou2) |laila 2 
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where we have introduced the compact notation Acting twice with the same sequence of operators gives 

IK)ar"K)„„) = \SZ 1 =s-mi)---\SZ n = s-m n ) 
gW.-aJS? = «)• (33) 



16(4s- 1) 



ii )1 r aia2 (1 <7 ai <7 a2 ) |2 ai 2 a2 ) 



(34) 



+ 



Er ai a 3 r a2a 3?7 aia2 .. . 
— ^ [(CT 0l + <T a3 )(Ta 2 ~ O ai <7 a3 - l\ |l 01 l oa 



4(4s-l) ^ 4s-r 0102 

a s I y/2s — 1 ^— \ r aiQ2 r aia3 ?7 Q2a3 . , M I o 1 1 

"Yj 2^ ~ 7TT + CT Q2 Cr Q3 - VaAVa 2 + <? a :i )\ |2 0l l 02 l 03| 



4(4s-l) ^ 8s-4 + T 02a , 



- y 

i6(4s - 1) ^ 8s - 2 + <7 ai o- a3 (r aia3 - r aia4 - r a2Q3 +r a 



-(1 - o- Ql er Q2 )(l 



where we have introduced the "non-coincident" symbol 
Vab = 1 - <W (35) 

The corresponding expressions for more successive ac- 
tions of these operators upon can also be obtained, 
but are too unwieldy to present here. 

Using such expressions, one may readily evaluate the 
terms, Hntlci}] in the diagonal effective Hamiltonian, 
Ea. (|2L)ll . For uq an even number, a convenient way to 
calculate the no-th order term is to consider the state 

|f) = i^^QHx [nQU^y 1 |-0 O ) (36) 

and then find the magnitude of this wavcfunction: 

n no [{ai}} = -<*|¥>. (37) 

Note that the square-root of 1Z in Eq. I)36|) is easily evalu- 
ated by just taking the square-root of Ea. (|31() . since it is 
diagonal in the basis of 3:1 configurations. Other terms 
can be obtained similarly. 

B. Restricting the Hilbert space to the 3:1 
manifold 

Calculating each such magnitude as defined in the pre- 
vious subsection leads to an explicit expression for the 
corresponding term in DPT. These expressions appear 
to be extremely complex and formidable functions of the 
Ising spin variables {c^}. In this subsection, we show 
that the projection of these functions to the 3:1 manifold 
of allowed {o~i\ configurations affords a tremendous sim- 
plification. In fact, we will demonstrate that all terms in 
DPT below 6th order can give only constant functions - 
i.e. no splitting - within the 3:1 states. At 6th order, 



the full functional dependence can be characterized by 
only 3 independent numbers which may be defined on 
plaquettes of the pyrochlore lattice. We show how these 
numbers can be extracted from the expressions obtained 
by the analysis of the previous subsection. 

1. Functional form of diagonal DPT terms 

From the analysis of the previous subsection, the n th 
order effective diagonal Hamiltonian in DPT must take 
the form of a multiple sum of n site indices a\---a n , 
where each site index is summed over all lattice sites. 
The summand is a function only of s and of the set of 
Sij where i, j must belong to the set of the site 
indices. The general form can be somewhat simplified 
by noting first that the dependence upon the r/y can 
be eliminating by rewriting them in terms of 5^ using 
Ea. (l35|) . and then eliminating all <5y by collapsing the 
sums containing these factors. Finally, we note that any 
factors of Ty in the denominators in these expressions 
can be moved to the numerator using the identity 

g(T ij )=g(0)+T ij [g(l)-g(0)} , (38) 

for any function g (which may also depend upon any 
other set of variables), since Ty = 0, 1. 

By these manipulations, one may write the effective 
Hamiltonian 7icff[{ci}] = J^n^niWi}] as 

H eS [{(Ti}} = (39) 

EE E II ^ kK,-^«J. 

n G„ ai—a-n\(ij)£G n J 

Here we have divided the effective Hamiltonian into terms 
involving n independent sites variables, a± . . . a n that are 




FIG. 1: Examples of contractible diagrams. 

summed over the lattice sites. A given order N in DPT 
contributes terms with n < N. For a given n, all pos- 
sible products of T aiaj can appear. The different such 
products are specified by G„, which may be considered 
as a "diagram" in the following fashion. Each G n can 
be represented by drawing n points, corresponding to 
i = 1 . . . n, and connecting some arbitrary set of pairs 
of these points by lines. For each (unordered) pair of 
points (ij) which is connected in G„, we include one fac- 
tor of T aia .. Since there are n(n — l)/2 pairs of points, 
and each pair may or may not be connected, there are 
2«(n-i)/2 d} s ti nc t diagrams G n . For example, in our con- 
ventions, r aia2 F a2a3 F a3a4 and r aia2 r a2a3 r a3a5 are rep- 
resented by different diagrams (see Fig.^l, which means 
that /g„(c q1 , . . . , cr a6 ) is not necessarily symmetric with 
respect to swapping a ai and a a& . We will refer to the 
number n as the order of the given term, even though it 
can come from a term of that order or higher in DPT. 

2. Contractible diagrams 

First we would like to show that any such term repre- 
sented by a diagram containing a point i with less than 
two connections to other points can be reduced to a term 
of one lower order. These diagrams are "contractible" 
(see Fig.|21for examples). We prove this by showing that 
the sum over cti can be carried out explicitly to obtain 
an expression of the same form of Eq. (|39|l in terms of the 
n— 1 remaining sum variables. There are two cases. Sup- 
pose in G n the point i in question has no lines connected 
to it. Taking i = n without loss of generality, we note 
that the sum on a n is unconstrained, i.e. it runs over all 
lattice sites. Thus we may write 

2 fc n (v ai ,---,°a n ) = J2J2f G " ^ ' ■ ■ ■ ' ( 4 °) 

a n t aEt 

= N t [3/ G „(cr Ql ,...,cr a7l _ 1 ,+) + / G „(cr ai ,...,cr Q „_ 1 ,-)] . 

The second line applies because on every tetrahedron 
there is the same set of four single-spin states. By in- 
serting Ea. (|40l) into Ea. (|39[) . one reduces the order of 
this term, as asserted above. 

Consider the second case, in which there is one connec- 
tion to the point i = n. We may suppose this connection 
is to the point j < n. The sum over a n is then con- 
strained only by the requirement that a n be a nearest- 
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1 1 1 

rr r£ 'Ci: 

(a) (b) (c) 

FIG. 2: Examples of contractible ((a) and (b)) and non- 
contractible (c) diagrams. 




(a) Site adjoining the (b) Two minority sites 

tetrahedra is the only connected by two 

minority site. parallel links. 




(c) Two minority sites 
connected by two links 
bending. 



FIG. 3: (Color online) The three possible configurations of 
minority sites (red arrow) on two adjacent tetrahedra, in the 
3:1 manifold of states. 



neighbor of aj. For fixed ay, this includes just 6 sites 
on the pyrochlore lattice. Moreover, the set of spins on 
these six sites is entirely determined by the spin at site 
Oj. In particular, if a aj = +1, the sum contains 4 terms 
with er Qri = +1 and 2 terms with o~ an = — 1; if a a . = — 1, 
the sum contains 6 spins with a an = +1. This can easily 
be understood from Fig. 01 Therefore the sum can again 
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•s m n 

FIG. 4: All n < 5 non-contractible diagrams. (a)The triangle 
diagram is the only possible such diagram for n = 3. (b) The 
square framed diagrams are all the possibilities for n = 4. 
(c) The pentagon-framed diagrams together with the three 
right-most diagrams comprise all the possibilities for n — 5. 

be carried out explicitly: 

y^r 0nCb ./g n (a ai) . . . ,<T a „) = 

On 

[4/ G „ (a Ql ,...,+) + 2f Gn (a ai ,...,-)] 
+^^6/ G „( ( r ai ,...,a a „_ 1 ,+). (41) 

Once again, Eq.l^lJ can be inserted into Eq.(j2Hl to re- 
duce the order by one. 

3. Non-contractible diagrams 

Since all contractible diagrams can be reduced using 
the above rules until they become either non-contractible 
or constant, we therefore need to consider only non- 
contractible diagrams. In these diagrams, each point in 
G n is connected to at least two other points. Let us 
first make a few general observations about these dia- 
grams. One can readily see that for these diagrams at 
order n < 5, all points must be connected, i.e. it is pos- 
sible to pass from one point to any other by a sequence 
of links. It is useful to consider the notion of a loop, or 
sequence of points, each connected to the next by a link, 
which visits no point twice and returns to the first point 
of the sequence. For n < 4, there is always at least one 
loop which includes all n points. For n = 5, all but three 
non-contractible diagrams contain a loop of length 5. The 
three remaining diagrams at n = 5 contain smaller loops 
(see part (c) of Fig. @J. All the non-contractible single 
loop diagrams for n < 5 are shown in Fig.0J For n = 6, 
there is one possible disconnected diagram, which con- 
tains two disjoint loops of length 3. Apart from this last 
diagram, all others are fully connected. 

Let us consider the physical pyrochlore sites which 
are summed over in a given term. They comprise a 
set S(G„) = {(aS 1) ,..., a [ 1) ),(af ) ,..., a l 2) ),...} of so- 
lutions, (fflj[ , ■ ■ ■ , On), to the conditions 

T aiaj = 1 for (ij) G G„. (42) 




(a) (b) <<0 



FIG. 5: (Color online) 3 example topologies of closed paths of 
five steps on the pyrochlore lattice. These are all the possible 
topologies of clusters corresponding to diagrams at order n = 
5 containing a loop of length 5. 

We will call these solutions "clusters" . In an infinite sys- 
tem, S is of course infinite because of translational sym- 
metry, but this is immaterial. A given term may then be 
written simply as 

f(cr ai ,...,a an ). (43) 

(oi,...,a„)eS(G„) 

We note that all the clusters for n < 5 are confined to 
one or two adjacent tetrahedra. This can be seen by 
considering the constraints imposed on clusters by the 
non-contractibility of the diagram. For instance, all but 
three diagrams at order n — 5 contain a loop of length 
5, and this allows only three topologies of clusters, which 
are illustrated in Fig. [5] The remaining three diagrams 
only allow clusters that are confined to two or less adja- 
cent tetrahedra. We will show more generally that any 
term containing only clusters confined to three or fewer 
adjacent tetrahedra is a constant. 

The set S can therefore be broken up into three com- 
ponents, comprising clusters which contain 1, 2, or only 
3 multiply-occupied tetrahedra, 

S(G„) = Si(G n ) + S 2 (G„) + S 3 (G„). (44) 

The sum in Eq. l|43|l can be carried out separately over 
these three sets. Let us consider first the sum over S±. 
The clusters in Si can be divided into subsets of those 
residing on a specific tetrahedron S* . 

An arbitrary permutation P of the 4 sites on tetrahe- 
dron t leaves the set Si(G n ) invariant. This is because 
each solution obeys Ea. H42(l . and T aiaj = rp( a .)p( aj ) for 
cii, cij G t (this is a set of permutations that leaves nearest 
neighbor pairs invariant). 

The contribution of all clusters on t to the term in 
question can only be a function of the 4 Ising variables 
of the 4 sites q = 1, 2, 3, 4 G t 

X] /G„(crai, ■ ■ -,VaJ = F(a 1 ,a 2 ,cr 3 ,a i ) . (45) 

(oi,...,a n )GS| 

Now we can use the fact that the spin configurations 
on one tetrahedron are always constrained to be of the 
3:1 form, i.e. they are a permutation P of the specific 
configuration + + + — : 

= (46) 
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with (ct?,o-^,<t§,ct2) = (+,+,+,-). Here q -> P(q) is a 
permutation of the 4 sites. The specific (cyclic) permu- 
tation P now encodes the spin state on this tetrahedron 



? \ a P(l) > a P(2) > °P(3) J °>(4) J 

X! fGn( a P(ai)> ■ ■ ' >°P(a„)) 
(oi,...,o„)eSf 

£ / G „«,..,<J 

(a 1 ,...,o„)eP-i(S t 1 ) 



(47) 



Since the set S\(G n ) is invariant under these per- 
mutations, we find from the last expression, that 
F{<J\, (72, <T3, CT4) is also invariant under the permutations. 
Hence this contribution is identical for all spin configu- 
rations, and is a constant within the 3:1 manifold. 

Let us next consider the clusters in S2. For each clus- 
ter, there are two neighboring tetrahedra t, if which each 
contain two or more sites . These tetrahedra share one 
specific site, which we call A. The pair of tetrahedra in 
question are determined by A (the tetrahedra t and t' 
are determined by the choice of A). For one such cluster, 
the sites ai with i = 1 . . . n may be partitioned into three 
groups: the site A, and those which are on t or if but are 
not A: 



t = {ai\ai e *,a» ^ A}, 
t' = {a 4 |a 4 e t',ai ^ A}. 



(48) 
(49) 



Similarly to Si, we can divide S2 into subsets S^ residing 
on tetrahedron pairs defined by the site A. We can then 
rewrite the sum by summing A over all lattice sites, and 
summing the set of sites a% . . . a n 



over S A 



E /G„(c r a 1 , ■ ■ ■ ,CT 0n ) 

(a 1 ,...,a 7l )dS2 

= E E fG n (<r ai ,---,<7a„ 
A (oi,...,a„)eS£ 



(50) 



We now observe that the set of solutions S2 is invariant 
under any permutation P t (P#) of the 3 sites in t (t'). 
Exactly as for S' this is because each solution in obeys 



Eq.g2Jl, and r 0<0j 



P t (o<)P t (a 3 ) 



for ai.aj S AUtUt' 



(and the same holds if Pt is replaced by Pt')- 
The sum 



E 

(ai,...,a„)GS2 



/G„(Cai) • ■ ■ , C a „) 



(51) 



can only be a function of the 7 Ising variables of the sites 
in A U t U t'. Due to the 3:1 constraint, if er^ = +, then 
the Ising variables a q for g £ t must be a permutation 

P t of a q 1] = (+ + -). If o a = -, then all the a q = +. 
Hence we may write 



(1 + tXA) (1) 

2 °> t (?) " 
(1+ CTA ) (1) 

~~ 2 a P t <(q) 



for g e t 



(+1) for get' 



(52) 



Using these expressions, and the fact that S A is invariant 



under these two permutations, 
found to depend only on oa- 
This leaves finally 



the sum in Eq. I|51|l is 



E =E/>^)' 



(53) 



where f{p a) is a complicated function obtained from the 
above manipulations - which however does not depend 
upon A itself. The sum is clearly then constant, as the 
number of + and - spins are fixed for the lattice. Thus 
all terms in S2 are also constants. 

Finally, consider S3. In these clusters there are three 
adjoining tetrahedra, and one may identify a "central" 
tetrahedron t which shares a site with each of the other 
two tetrahedra if, t" . Here one may divide the sum vari- 
ables into five groups: two corresponding to the site 
shared by t, t' and the site shared by t, t" , and three oth- 
ers corresponding to the sites on i, If, t" but not shared. 
One can again sum over the unshared sites on t' and 
t", and obtain an expression for the cluster sum which 
involves sites only on t. By manipulations of the type 
used to analyze Si, one finds that this remaining single- 
tetrahedron sum must also be constant. 

We conclude that any term for which the correspond- 
ing clusters are confined to three or fewer adjacent tetra- 
hedra must be constant. Therefore all terms up to and 
including 5th order are constant. At sixth order, amongst 
the non-contractible diagrams there are a few exceptions. 
First, there is one disconnected diagram containing two 
loops of length three. In this term, the sum over vari- 
ables in the first and second groups is independent, and 
therefore each can be carried out separately as for a third 
order term. This gives immediately a constant contribu- 
tion. The remaining diagrams are connected. All but one 
of these diagrams contains a loop of length 5 or less (pos- 
sibly in addition to other larger loops). Such terms are 
confined to three or fewer tetrahedra, and are constant 
by the above arguments. What remains is the single dia- 
gram consisting of only a single loop of length six, shown 
in Fig© 

This "large loop" diagram is thus the sole non-trivial 
contribution. It can be written in the form 



n^[{<7 l }} = 



ai...a 6 



where we identify ay — a\. To analyze each term (given 
a particular set a\ . . . as), we employ a trick: multiplying 
it by a carefully chosen representation of the identity 



1 



n c 



(55) 



with r] a ]y = 1 — 8 a b- Here the product is over distinct 
pairs i,j which are not connected in the loop diagram. 
We multiply the loop term by this expression and expand 
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FIG. 6: The only diagram at order n = 6 giving a non- 
constant diagonal contribution in degenerate perturbation 
theory. 



the product fully. All but one term involves at least one 
Kronecker (5-function. In each of these summand terms, 
at least one sum can be collapsed, leading to a lower- 
order term, which we have already shown is necessarily a 
constant. The remaining non- vanishing part is the origi- 
nal summand term multiplied by the product, 



(56) 



This factor is non-zero if and only if all n = 6 sites Oj are 
distinct. Thus the sites a, must comprise a closed walk on 
the lattice in which each site on the walk is visited only 
once. On the pyrochlore lattice, this is exactly the set of 
hexagonal plaquettes. A specific plaquette on the lattice 
containing sites i\ . . .is in sequence around the plaquette 
appears 12 times in the sum in Ea. (|54fl . with a\ . . . a§ 
taking the six cyclic permutations of i\ . . . i% and the six 
cyclic permutations of these sites in reverse order. Hence 
the non-constant contribution to the diagonal energy at 
6th order in DPT can be written: 



He, 



V 



&p(ii 



(57) 



where i±, . . . , i§ are the six sites moving clockwise around 
plaquette V, and 

£p(fiii • ■ ■ i fi 6 ) = 
6 

[h(<7i k ,- ■■ ,^ k+5 ) + f L (cr lk+5 ,.. .,0-iJ] , (58) 

k=l 

where i k+6 = i k . 



DPT calculation gives a specific energy (proportional to 



Type 


Configuration 


Fraction of 
minority spins 










1 


mm 


i 

2 


2 


ITTTTT 


1 

6 


3 


mm 


1 
3 


4 


mm 


1 
3 



TABLE I: The different plaquette types, with the fraction of 
minority sites in each one. 



There is some freedom in the choice of these 5 energies. 
That is, certain changes of the plaquette energies leave 
the differences of total energy amongst distinct 3:1 states 
unchanged. One obvious such "gauge" change is a global 
shift of all 5 energies by the same amount. Another less 
obvious constraint comes directly from the 3:1 rule. If 
one denotes the fraction of plaquettes in the lattice in 
configuration a by x a , the total fraction of minority sites 
must always be 1/4. Each plaquette configuration has 
a fixed fraction of minority sites M a , given in Table [I] 
Thus 



M a x a . 



(59) 



a=0 



The energy per plaquette is then 



^6 ^ ^a^a 



(60) 



a=0 



Using H59|l. one sees that a shift A£ a = cM a , with ar- 
bitrary constant c shifts the energy by a constant. The 
obvious global energy shift remarked on above derives 
similarly from the normalization condition ^ a x a = 1. 
Using these two constraints, we see there are only 3 in- 
dependent plaquette fractions. We (arbitrarily) choose 
to keep £1,2,4 as our independent variables. Substituting 
the solutions for the other fractions (0:0,3) into l|6"U|) . we 
find 



C. Results 

We have carried out the calculations detailed in the 
previous subsections. Specifically, by explicitly con- 
structing I*) in Eq.JSSJl, we obtained H 6 [{oi}] in Eq. (|3Tjl . 
From this, we extracted the function Jl in Ea. H54|) 
and thereby determined the plaquette energies E-p using 
Ea. (|58fl . Using the 3:1 constraint, there are 5 configura- 
tions possible on any plaquette, which we denote "type 
0" to "type 4" . These are enumerated in Table |]] The 



Ha = x x Vi +x 2 V 2 + x 4 V 4 , (61) 
with the 3 "gauge invariant" physical energy parameters 

Vi = i {£ + 2£ 1 - 3£ 3 ) , 

V 2 = ^(-£o + 2£ 2 -£ 3 ), 

Vi = (£i-£ 3 ). (62) 
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Our DPT results are: 



j Q 6 3s4 ( 983 04s5 ~ 139648s 4 + 79136s 3 - 22040s 2 + 3006s - 165) 



32(2s - l)(4s - l) 5 (8s - 3) 2 (12s - 5) 



V 2 =J z a { 



Va =J 7 a b 



s 3 (256s 3 - 51s + 9) 
32(4s- l) 3 (8s-3) 2 ' 
s 4 (272s 2 - 136s + 15) 
16(4s- l) 5 (8s- 3) 2 



(63) 



r 



We have made several checks on the above calculation. 
First, we have carried out a more direct scheme which 
sums the terms in DPT in a completely different manner 
from the methods described in this section. We leave 
the vast details of this calculation to Appendix [51 The 
results of this alternative method agree perfectly with 
those quoted above. Second, in the following section we 
will compare the s — > oo limit of the above result with 
the result of a large s calculation for the XXZ model. 
The large s limit of the energies we find in DPT becomes 



lim — =0, 

s — >oo s 



lim 



V 2 J z a e 



s~ >oo s 512 
Va 

lim — =0 . 

S — >OQ S 



(64) 



We shall see that this result indeed agrees exactly with 
the corresponding limit of the large s expansion. 



Each one of the DPT processes contributing to the 
off-diagonal term consists of 2s spin transfer operations 
along each one of 3 links of a hexagonal plaquette of the 
pyrochlore lattice (see Fig. 01 , acting in some particular 
order. 

We can calculate K by adding the contributions from 
all the DPT processes occurring on a single plaquette, 
starting in the state | ITITIT) an d ending in the state 

I nun). 

In every such process 3 spins go from an initial state 
of +s to — s, and 3 start with — s and end up as +s. 
The spins change via ladder operators 5~ , and therefore 
we get "angular momentum factors" from the action of 
these operators. The same set of operators S 1 * act in 
every process, and so these factors are always the same. 
For the S + operators taking a single site from — s to +s 
we find 



j j y/s(s + 1) - m(m + 1) = (2s)! 



(66) 



D. Off diagonal term 



and for the S operators taking a single site from +s to 
— s we find 



In this section we describe how the lowest order off 
diagonal term in the DPT effective Hamiltonian is calcu- 
lated. As explained in Section fill A II this term appears 
at order 0{a 6s ). 

The lowest order off diagonal term acts only on a 
hexagonal plaquette in the flippable configuration (type 1 
plaquette, as in Table fj}). It changes the plaquette con- 
figuration from one flippable configuration to the other 
flippable configuration. Therefore the off diagonal term 
has the following general form 

Hoff diagonal = (-1) 6S+1 a 6s J z K ^ ^' ^ ^ + ^ 

P 

=(-i) es+i a 6s j z i<J2(\ mmxnnn i + h.c) , 

(65) 

where we denote the two flippable configurations of the 

plaquette by ! •-*), ("-• I for the sake of brevity. Note that we 
can change the (— l) 6s+1 factor into a (—1) by a unitary 
transformation similar to that employed in Ref. 0- We 
shall now describe how the coefficient K is calculated. 



J] y/s(s + 1) - m{m - 1) = (2s)! . (67) 

m— — s+1 

In total from all the ladder operators, we find a common 
factor ((2s)!) 6 . From the 6 spin transfer operators we 
have another common factor of l/2 6s . 

All that remains to be calculated for a single DPT pro- 
cess is the product of resolvents of each stage in the spin 
transfer process. First let us classify the different pro- 
cesses on a single plaquette. We can choose one of two 
sets of three links on which spin transfer will occur (one 
such choice is shown in Fig. EJ) ■ The contribution from 
each one of these two cases is identical, so we shall calcu- 
late the contributions for one set of three links, and mul- 
tiply the final result by 2. The processes we will sum over 
only differ by the order in which the spin transfer opera- 
tors act on the 3 predetermined links. We call the three 
links A,B, and C, and then each process is described by 
a string of 6s letters qi . . . qe s which contain 2s instances 
of each one of the three letters A,B,C. For example, a 
possible string for s = 1 is AABBCC. From this classi- 
fication, it is evident that in total there are ( 2 s)!(2s)!(2s)! 
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then rewrite the resolvent as 



FIG. 7: (Color online) Off diagonal process on a single pla- 
quette. The (red) circles denote minority sites. 



different processes. At this point we can write a formal 
expression for the coefficient K 



6s- 1 



K 



— 2^ 11 n dUn}) , 



2 6s 



(68) 



{g„} i=l 



where Tte({Qn}) denotes the resolvent at step I of the 
DPT process described by the string {(?„}. 

Now we turn to formulating the resolvent in a con- 
venient manner that will facilitate the summation over 
all processes. Starting from Ea. (|31[) . in this case the set 
F consists only of the 6 sites surrounding the hexagonal 
plaquette. We shall denote these 6 sites 1 through 6, as 
in Fig. so that A denotes the link (1,2), B denotes 
the link (3,4), and C denotes the link (5,6). Since the 6 
sites have alternating initial states ±s, any pair of nearest 
neighbor sites has a^aj = — 1. We can therefore rewrite 
the inverse resolvent operator as 



■R- 1 = — J- 



6 

£ 



(SJ-s) (5| +1 -s)-2J, 



6 

»£ 

3=1 



(69) 

where the indices are defined modulo 6, so that Sf +1 = 
Sf. From this point on, all index arithmetic is defined 
modulo 6 as well, for ease of presentation. 

To further simplify the resolvent, we introduce 
iia{£, {q n }) as the number of times the link A has had 
spin transfer occur on it up to stage £ in the process 
desribed by the string {q n }. The same numbers can be 
also introduced for B and C. Then, by definition, the 
total number of spin transfer operations is Ua(£, {?n}) + 
n B{£,{q n }) + nc(£,{q n }) = £■ In what follows we will 
show that the resolvent can be described only by these 3 
numbers. To see this, notice first that, regardless of the 
order of spin transfer operations, a spin transfer operator 
on the link (j, changes (S*| — s) — ► (S| — s — l) and 
(Sj+i — s) — > (Sj +1 — s — l) . Note also that, in the ini- 
tial state, all (<S| +1 — s) = 0. Thus, at every stage of any 
process, (Sf — s) — (Sf — s) = — tia(£ , {q n }) Similarly 
(Sf-s) = (Si-s) = -n B (l,{q n }), and (Sf - a) = 
(iSf — s) = —nc{i, {q n })- Using these variables, one can 



72*({?„}) = 4J z sl 

_ 3 JL 
2 



(tia + ub) + {nB + nc) 2 + [nc + nA) 2 



(70) 



where we have suppressed the explicit dependence of the 
n A,B,c numbers on t, {q n } for clarity, ft is more conve- 
nient to derive a recursion relation for the resolvent at 
stage t 

ne+i({q n }) = Ki({q n })+Jz (4s - 1 - £ - n qe+1 (£, {q n })) 

(71) 

The initial condition for this recursive series is TZ 1 = 
0. Using Eq.J7TJ, we can calculate the product 
Ilfci 1 T^-eiiln}) for a given process. For every process, 
we need to keep track of only the 3 numbers ua.b,c m 
the various steps of the the process. 

We have calculated the coefficient K explicitly for a 
number of interesting values of s. The results are sum- 
marized in Table ITU 



s 


K 


1 


3 


2 


2 


1 


0.166 


3 
2 


0.00839536 


2 


0.000304464 


■ r > 

2 


9.1747 x 10" 6 



TABLE II: Values K of the coefficient for the lowest order 
off-diagonal term, for various values of s. 



IV. LARGE s EXPANSION 

A large s analysis has recently been employed in Ref . 26 
to explore the magnetic order for the general spin s 
Heisenberg AFM on the pyrochlore lattice. Restricting 
the Hilbert space to collinear spin configurations, the au- 
thors of Ref. derived an effective Hamiltonian out of 
the harmonic spin wave energy contribution, to order 
0(s). The effective Hamiltonian prefers spin products 
around hexagonal plaquettes to be —s 6 in the zero mag- 
netic field, and +s 6 in the half-polarized plateau region. 
Following a terminology inspired by Ising gauge theory, 
these are denoted u ir flux" configurations and "zero flux" 
configurations, respectively. In order to compare this ap- 
proach, which is justified in the large s limit, with the 
DPT analysis of Section IfTTI we have repeated the same 
type of effective Hamiltonian calculation for the the XXZ 
model. Our derivation follows closely that of Ref. 12ft 

The large s expansion consists of expressing the spin 
degrees of freedom in terms of Holstein-Primakoff bosons 
and expanding in decreasing powers of s. The lowest 
order term in the large s expansion is of order s 2 , and 
corresponds to the classical spin version of the quantum 
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FIG. 8: (Color online) Magnetization (in units of s) of a single 
tetrahedron of classical spins with an anisotropic XXZ inter- 
action, parametrized by a. For any a < 1 a half polarization 
plateau exists. 



XXZ Hamiltonian 
«d =J Z [« (Si • Sj) + (1 - a) (Si ■ z) (Sj ■ z)] 



2J z h J2 S ! > 



(72) 



where as before a — J±/ J z . 

In order to analyze the ground state of this anisotropic 
classical model <|72[l . we first calculate the minimum en- 
ergy configuration for a single tetrahedron. For the single 
tetrahedron we obtain the magnetization curve shown in 
Fig. [5] We find that for a < 1 a half polarization plateau 
opens up, and the plateau becomes wider as a decreases 
from 1. In this plateau, the classical spins on the single 
tetrahedron being analyzed are in a collinear configura- 
tion, with three Sj = sz and one Sj — —sz spins. This 
is just the classical analog of the 3:1 configuration on a 
single tetrahedron found in Section^ A 3:1 spin config- 
uration can be realized on each and every tetrahedron of 
the lattice simultaneously. We therefore conclude that in 
the range of magnetic fields where the single tetrahedron 
is in a half polarized state, the ground state of the many 
body system (I72[) is the manifold of 3:1 configurations. 
This means that the plateaus in the classical XXZ model 
on the complete pyrochlore lattice are at least as wide as 
in FigEl 

In the following we will discuss only this half magneti- 
zation plateau. We then assume the collinear 3:1 states, 
which allows us to describe the magnetic configuration in 
terms of the same Ising variables crj = ±1 as in Sectionlrfl 

As in Ref . I2t3 we use the unitary transformation l|17|) so 
that we can define the Holstein-Primakoff bosons, which 
amounts to replacing the rotated spin operators as follows 



S| =» - 



S+ =v /2s - rhj bj 



'2s bi 



(73) 



where bj are canonical bosonic operators, and rhj = b^bj 
is the boson number operator. We plug these into the 
Hamiltonian and keep only the quadratic terms in 
the bosonic operators. 

Since the spin configurations are now restricted to the 
3:1 manifold, the magnetic field term is the same for ev- 
ery 3:1 configuration as the magnetization is constant on 
the plateau. In terms of the unrotated spin variables S* , 
this amounts to J^j ^j = i^ where N is the number of 
sites in the pyrochlore lattice. Varying the magnetic field 
in the plateau region causes an overall shift in the spin 
wave energies of all the 3:1 states, and thus will not alter 
the energy differences between different 3:1 states. Simi- 
larly the Ising variables have a sum of J2j a j = 2^> anc ^ 
we can use these two identities to derive ■ Ojihj = 0, 
which is useful in simplifying other terms. Therefore, we 
can ignore the magnetic field term, since we are search- 
ing for an effective Hamiltonian splitting the energies of 
different 3:1 states. The effect of the magnetic field is to 
determine the energy gap for spin wave excitations. The 
vanishing of the spin wave gap signifies an instability of 
the 3:1 manifold, corresponding to the edges of the half 
polarization plateau. 

From Eas. (|B3|l . 1)23(1 . the resulting harmonic spin wave 
term reads 



o-3:l j f* r 

rl harm J z „ b / J 1 



1 + OjCT 



■)( 



b]h + h.c. 



1 — (TjCTj 



b-jbi + h.c. 



J z 2s^TOj. (74) 



Following the derivation Ref. l2(:l the zero point energy 
of this harmonic term for a given 3:1 spin configuration 
(described by {(Jj}f =l ) is 



N 



(75) 



k=l 



where Afe are the solutions of the eigenvalue equation 



(*f* + f) + (!*f) : 



v . (76) 



In the right hand side f denotes the same N x N connec- 
tivity matrix introduced in Section^ and a is a diagonal 
N x N matrix with <jj as its diagonal elements. Without 
specifying the 3:1 configuration, we can write an expres- 
sion for the harmonic energy in terms of Uj 



hari: 



J z sTr 



1 + £ (ite + f) + f(,f)- 



(77) 

One can calculate the spin wave energies by assuming 
a particular spin configuration and computing the trace 
exactly. However, as in Ref. 12(1 if one does not know 
which candidate spin configurations to consider, one can 
derive an effective Hamiltonian to determine which spin 
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configuration gives the lowest harmonic energy, and find 
a favorable spin configuration. 

The square root in Q77JI can be expanded in powers of 
matrix operators. We first observe that a only appears as 
a multiplier of the matrix T. Therefore, an expansion in 
powers of matrix operators is equivalent to expansion in 
the parameter a. In the present context, this expansion 
is justified due to the easy axis anisotropy a < 1. 

The terms in the expansion can be organized as a sum 
of traces over products of T matrices and Ising variables 
tjj . The order of a for each term also specifies the number 
of connectivity matrices T appearing in that term. 

Due to the trace operation, the product of T matrices 
represents closed loops on the lattice. The Ising variables 
appearing in each such term can only involve the sites on 
the loops defined by the product of T matrices. Using the 
results of Section AlI B 21 which discuss functions of Ising 
variables and T matrices precisely of the form appearing 
in this expansion, it is evident that all terms involving 
less than six T matrices will result in constants, which will 
not split energies of the 3:1 states. As in Section the 
lowest order term in the expansion in a causing energy 
splitting in the 3:1 manifold involves loops around hexag- 
onal plaquettes of the pyrochlore lattice. For simplicity, 
we consider only these terms, and ignore any higher order 
term in the expansion in a. After extensive simplifica- 
tion, the 6-th order term reads 



a 



" ( 2 ) 512 

+ 7Tr (a.T 3 .c7.T 3 ) - 6Tr (cr.T.a.T.cr.T.a.r 3 ) 



14Tr (a.r.a.T 5 ) + 14Tr (a.T 2 .a.T 4 



3Tr(a.r 2 .cr.r.cr.r 2 .cr.r) 
6Tr(a.r 2 .cr.r 2 .cr.r.cr.r) 

Tr(CT.r.cr.r.ff.r.CT.r.CT.r.cr.r)] + o{a s ) 



(78) 



From the expression one extracts only those terms 
corresponding to loops around hexagonal plaquettes. 
Eci. (|7Sj) takes the form of the function in Ea. (|3^|l . 
with n = 6 and the "loop" diagram G n — 
{(12), (23), (34), (45), (56), (61)}. The corresponding 
function /(c ai , . . . , o~ a6 ) reads 

+ 14 <7 ai cr a3 + 7 a ai a ai 

— 6 <Ta 1 <Ta 2 &a 3 &a i — 3 <7 ai CTa 3 CTa 4 (J a & — 6 (T 0l a a3 <Ta 5 O a & 

(79) 

The effective Hamiltonian therefore describes all possi- 
ble spin interactions on the hexagonal plaquette of the 
pyrochlore lattice - 2,4 and 6 spin interactions. It is far 
more convenient to express this complicated Hamiltonian 
in terms of energies of plaquette configurations, in the 
same way we formulated the results of the DPT in Sec- 
tion [n] as Hharm = J2p^p (using the same 5 plaquettes 

in Table 0. 

As in Section Till CI there are only 3 independent pla- 
quette configuration energies, V\^,a which to 0(a 6 ) are 



Overlap 




FIG. 9: This figure shows the regions of parameter space 
where the DPT and large S expansions are justified, and their 
region of overlapping validity. 



Pi =o, 

J z a 6 s 



V 2 =- 



512 



(80) 



Va =0 



Comparing H8U|) with (|64|l , we find complete agreement 
between the DPT of Section UTI and the large s expansion 
of this section, in the limit of both a — > and s — > oo, 
where both approaches are justified (see Fig. |!5J). This 
serves an excellent check on the correctness as well as va- 
lidity of our calculations, in the parameter regime where 
the approximations overlap. 



V. LOW ENERGY STATES OF THE 
EFFECTIVE HAMILTONIAN 

A. Strict easy axis limit for s > 3/2 

In this subsection, we consider the a — J±/J z <C 1 
limit, for which the lowest-order non- vanishing terms in 
the effective Hamiltonian are dominant. For any s > 3/2, 
this is just the sixth order diagonal contribution. 



1. large-s case 

We first consider the large s limit. As is clear from 
Eq. I|8U|I , at order s only the type 2 plaquette suffers from 
a positive energy correction, while Pi and V4 are only 
nonzero at order s° or lower. Hence the type 2 plaquette 
is strongly disfavored for large s. This in combination 
with the 3:1 constraint allows us to restrict ourselves to 
the "0-flux manifold" in the large s region (See Sect ion llVl 
for the definition of a 0-flux manifold.) . 

To see this, let us first introduce a "cell" comprised by 
4 link-sharing plaquettes. Choose 4 hexagonal plaquettes 
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such that any pair of two plaquettes out of these 4 pla- 
quettes always share a link. Then these four plaquettes 
form a single polyhedron (a truncated tetrahedron), with 
4 hexagonal faces, and 4 triangular faces (see Fig. I10[) . 
We will refer to this polyhedron as a cell. In the py- 
rochlore lattice, one may distinguish two kinds of cells - 
when one completes the tetrahedra enclosing a cell, we 
can identify up-headed and down-headed cells, according 
to the direction at which the tetrahedra are pointing (see 
Fig. ^1 for examples of both kinds ). Each up/down- 
headed cell shares its faces (hexagonal plaquettes) with 
4 nearest neighboring down/up- headed cells. Thus, cen- 
ters of cells constitute a diamond lattice, where those of 
up-headed cells take part of one FCC lattice and those 
of down-headed cells form the other FCC lattice. It suf- 
fices to determine the spin configuration on only the up- 
headed (down-headed) cells in order to specify the spin 
configuration on all sites of the lattice. 

Observing the local constraint, one can readily enu- 
merate the various minority spin configurations of a cell. 
In Table. 11111 all possible cell configurations allowed in 
the 3:1 manifold are listed. Each cell type is described 
by the configurations of its 4 hexagonal plaquettes. 

To see that the ground state manifold in the large-s 
region is composed only of type 0, 3 and 4 plaquettes 
(i.e. 0-flux states), notice from Tabic. ITTT1 that any cell 
type which contains a type 1 plaquette always contains 
at least one type 2 hexagonal plaquette. This implies 

< xi < x 2 (81) 

where the plaquette type fractions, as introduced 

in Section IIII CI Disallowing the type 2 plaquette in- 
evitably leads to excluding the type 1 plaquette, and 
therefore, the positive Vi in leading order of s _1 expan- 
sion allows us to conclude that the classical ground state 
spin configurations in the large s limit consist of only the 
0-flux states. 

This ground state, however, is massively degenerate. 
Higher order quantum corrections in s _1 can select a par- 
ticular classical state out of this 0-flux manifold. To see 
this, let us expand the plaquette energies in 

Vi 17 ( -2\ 
J z a 6 ~ 65536s + ^ ' 

Notice first that the 0(1) negative energy correction to 
V\ plays no role in lifting the degeneracy of the 0-flux 
manifold, since this manifold does not contain any type 1 
hexagonal plaquettes. Thus, provided that V% dominates 
the other two, the most relevant correction in the large s 
limit is V4, which always disfavors the type 4 hexagonal 
plaquette, since it is positive. 

Since the type 4 plaquette is disfavored, observing the 
0-flux condition on all plaquettes, we have only to mini- 
mize X4 to obtain the ground state in the large s region. 
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TABLE III: The various cell configurations are described by 
the number of each plaquette type included in the plaquettes 
comprising a cell. Cell types are indicated by italics and pla- 
quette type by sans serif. In the zero-flux manifold, only the 
type 9, 10, 11 and 12 cells are allowed, since they do not con- 
tain type 1 and type 2 hexagonal plaquettes. Furthermore, a 
cell must contain a type 2 plaquette whenever it contains a 
type 1 plaquette. This is quantified by < x\ < X2- 




FIG. 10: (Color online) A cell is composed of four link- 
sharing hexagonal plaquettes (The polyhedron bounded by 
thick (blue) lines). The cell on the left is up-headed, and the 
one on the right is a down-headed cell. 

However, in the 3:1 manifold, the 0-flux condition be- 
comes so strong that X4 is in fact bounded by J5 from 
below (X4 > Jr). To see this, notice first that only 
the type 9, 10, 11 and type 12 cells drawn in Figllll 
are allowed in the 0-flux manifold. Next, we denote by 
2/9,10,11,12 the fraction of cell types 9 ... 12 in the entire 
pyrochlore lattice (we use these instead of plaquette type 
fractions (x^, X4, xo) for later convenience). In the 0-flux 
manifold, only these cell types may occur, and therefore 
Y^jLgUj = 1- Together with the "global" 3:1 constraint, 
i.e. Eq. (|5^|l one finds 

3yi2 = 2/9, (83) 

or alternatively, 

2/12 = j(l - 2/10 - 2/11), ■ (84) 

An important step to identify the lower bound on £4 
is to note that packing these four cell types into a py- 
rochlore lattice is highly constrained by the local 3:1 rule 
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type 9 cell 



type 10 cell 




type 11 cell 



type 12 cell 



FIG. 11: (Color online) The four cell types allowed in the 
zero-flux manifold. Minority sites are specified by the (red) 
circles. 



imposed on each tetrahedron. For example, a type 12 
cell can only have cell types 10,11 and 12 as neighboring 
cells. Each type 10 and 11 cell can neighbor at most one 
type 12 cell, as they both have only one type plaquette, 
and the type 12 cell consists only of type plaquettes. 
One can also show that a type 12 cell can have at most 
one neighboring type 12 cell, the remaining neighboring 
cells must be of type 10 or 11. These observations are 
already sufficient to conclude that 



3yi2 < yio + 2/11- 



(85) 



Now using Eq. JH3J, we obtain the lower bound on the 
fraction of type 10 and type 11 cells: j/io + 2/n > f • 
Since these two types of cell are the only cells allowed in 
the 0-flux manifold which have type 4 hexagonal plaque- 
ttes, this lower bound immediately gives us that for the 
fraction of type 4 hexagonal plaquette: 

13 113 
%i = -jVio + ^2/n > j2/w + j2/H > ^g- (86) 

From the derivation above, one can easily see that the 
equal sign is realized if and only if yu = 0. Excluding 
type 11 cell configurations, one can show that a type 12 
cell always neighbors three type 10 cells, and one type 
12 cell. As a consequence, the condition 3yi2 = 2/10 is 
satisfied only when any type 10 cell has a type 12 cell 
as its nearest neighboring cell, through its single type 
plaquette. Since the fraction of type 9 cells is uniquely 
determined by the fraction of type 12 cells (|83l) . and we 
have already excluded any type 11 cells from a state sat- 
urating the bound on if a state saturating this bound 
exists, since it minimizes the fraction x± it must be con- 
figuration with a maximum number of type 12 cells on 
the lattice. Without any type 12 cells, we cannot have 




FIG. 12: (Color online) The trigonab state. The spin con- 
figuration of a planar layer of tetrahedra is shown. Triangles 
with lines connected at their centers represent up pointing 
tetrahedra, while the other triangles represent down pointing 
tetrahedra. Minority sites are denoted by (red) circles. Two 
dashed triangles denote up pointing tetrahedra in the planar 
layer of tetrahedra immediately above the one depicted in this 
figure. These two tetrahedra are used to show the primitive 
vectors for the pyrochlore lattice ai,2,3 and for the magnetic 
unit cell of the trigonab state Ei,2,3- The seven up pointing 
tetrahedra included in one valid choice of a magnetic unit cell 
for the trigonab state are marked by (blue) letters indicating 
one of 4 3:1 configurations for a single tetrahedron. The type 
plaquettes residing between pairs of adjacent type 12 cells 
are marked by large (blue) circles. 



any type 9 cells either, and are limited to type 10 and 11 
cells. A state comprised only of type 10 and type 11 cell, 
has fraction of type 4 plaquettes that is always greater 
than |. 

In what follows, we will show that this lower bound for 
X4 is uniquely (up to a finite degeneracy) realized by the 
periodic minority spin configuration depicted in Fig. 1121 
This collinear magnetic ordered state, which we shall re- 
fer to as the "trigonab" state, contains 7 pyrochlore unit 
cells. It has a magnetic unit cell with primitive vectors 
Ei = 2ai — a3, E2 = 2a2 — ai, and E3 = 2a3 — a2, where 
ai.2,3 are the primitive unit vectors of the pyrochlore lat- 
tice (FCC lattice vectors a! = #-(0, 1, 1) and cyclic per- 
mutations). From the unit cell vectors, we can find the 
volume of the magnetic unit cell 



(Ei x E 2 ) • E 3 = 7 (ai x a 2 ) • a 3 



(87) 



These 3 primitive vectors are of equal length, and are not 
mutually perpendicular. Therefore, the magnetic Bravais 
lattice is in the trigonal crystal system - whence the name 
trigonal state. 

From the planar view in Fig. ^| it is clear that this 
magnetic state has a three-fold rotation symmetry about 
the ai+a2+a3 = a(l, 1, 1) axis perpendicular to the page. 
In the direction of the pyrochlore lattice directions, there 
is a periodicity of 7, giving rise to a seven fold degeneracy 
due to FCC lattice translations alone. The trigonab state 
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FIG. 13: (Color online) The broken reflection symmetry. 
Starting from a type 12 (up-headed) cell, and drawing the 
10 tetrahedra surrounding it, we first choose the place of the 
nearest neighboring type 12 cell (down-headed) - the two type 
12 cells share the hexagonal plaquette marked by thick dashed 
(blue) lines. With this choice, the minority sites on 7 tetra- 
hedra are automatically determined (marked by open (red) 
circles). However, minority sites for the other 3 tetrahedra 
(solid (green) circles) are not fully determined and we still 
have a "mirror" degree of freedom. For convenience, we also 
draw the primitive vectors of pyrochlore lattice, in accordance 
with those defined in Fig. 1121 



breaks a reflection symmetry about a plane perpendicular 
to the Kagome plane, parallel to ai and passing through 
the point where the three vectors Ei,2,3 originate in the 
figure (see Fig. 1131 for a another view of this symmetry 
operation). Together with the 4- fold choice of the set of 
Kagome planes, it is evident that the degeneracy of this 
magnetic state is 4 x 7 x 2 = 56 

As is clear from Fig.^J the spin configuration satisfies 
both the local zero flux condition and the local 3:1 con- 
straint. To sec that this trigonal state saturates the lower 
bound for x±, one has only to count the fraction of type 
9, 10, 11 and 12 cells: yg : y%o : yn : yx2 = 3 : 3 : : 1. 
From Table IIIII we find the trigonal state realizes the 
lower bound X4 = We conclude the trigonal state is 
at least one of the ground states in large-s region. 

As argued above, any state saturating the bound must 
have every type plaquette connecting between a type 
12 cell and a type 10 cell. Starting with a type 12 cell, 
and using this rule together with the 3:1 constraint and 
the zero flux condition, suffice to uniquely construct the 
trigonaly state, up to the finite degeneracy described in 
above. Starting from the initial type 12 cell, the plaque- 
tte connecting this cell to another type 12 cell defines the 
Kagome plane in Fig. ^] Next, pick one of the 2 mirror 
image choices in Fig.^|of the type 10 cell configurations 
neighboring the first type 12 cell. From this point on, 
the three rules mentioned above uniquely determine the 
rest of the magnetic configuration in the entire lattice. 

Finally, the energy per plaquette of the trigonaly state 

is 

1 3 

■rr-Btrigoiial 7 = 7^ V 4 ■ ( 88 ) 



2. Spin s>2 

We expect that the trigonaly state described above is 
the ground state for sufficiently large s. In the following, 
we shall argue that this is indeed the case for s > 2. For 
s = 5/2, 2, 3/2, 1, the energy parameters in the effective 
Hamiltonian are given in Table ITVl 

For all the cases in Table IIVI V\ is the largest and 
most negative energy. This would suggest that the low- 
est energy 3:1 state is one with a maximum number of 
plaquettes of type 1. However, the geometry of the lattice 
as well as the 3:1 constraint pose stringent restrictions. 
By enumerating all possible types of cells in Table IIIII 
one finds that every type 1 plaquette must be accompa- 
nied by at least one type 2 plaquette on the same cell. 
The configuration of the entire lattice can be determined 
by considering only up-headed cells, and therefore, the 
existence of M type 1 plaquettes demands that at least 
M type 2 plaquettes are present as well. We deduce the 
following inequality x\ < x%, in any 3:1 configuration. 
From this we see that the energy of a type 1 plaquette is 
offset by the energy cost of a type 2 plaquette which is 
the highest energy cost for all the s values in Table Hvl 
Therefore, the number of type 1 plaquettes is not nec- 
essarily maximized in the ground state even with small 
s. 

One observes that the magnitude of the energy V\ is 
comparable to V2 already at s = 5/2, and this trend 
continues to higher s -V2 becomes more dominant. Given 
the restriction x\ < X2, and the large energy cost of 
type 2 plaquettes, the analytic arguments in the above 
subsection suggest the trigonaly state may be the lowest 
energy state for s > 5/2. The case s — 2 is close to the 
boundary for a change in behavior. 

In order to search for other candidate ground states, 
we have enumerated all 3:1 states on a variety of peri- 
odic finite clusters, and determined the exact lowest en- 
ergy state for each one, for s = 1,3/2, 2, 5/2, • • • , 6. For 
s > 2 we find no states with lower energy than that of the 
trigonaly state. This strongly suggests that the trigonaly 
state is the ground state for all s > 2, though of course 
this limited numerical investigation does not constitute 
a proof that this is the case. Moreover, states with large 
numbers of type 1 plaquettes, are among the highest en- 
ergy states we have found, which does give credence to 
our assessment that when the V% and V\ are comparable 
energy scales (with opposite sign) , because of the condi- 
tion x\ < X2 the energy V2 is still dominant. One can 
conclude that if there is a state with lower energy for 
s > 2, it must have a large unit cell which is incompati- 
ble with all the clusters considered in Table Ivl 

3. Spin s = 3/2 

Spin s = 3/2 is the smallest spin value for which in 
the extreme easy-axis limit a <C 1 the off-diagonal term 
in the effective Hamiltonian may be ignored. The cor- 
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TABLE IV: Energies Vi,2,4 of the plaquette configurations 
type 1,2,4 for s = 2, §, 1 
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TABLE V: 3:1 configurations on periodic clusters. Energy is 
given in units of J z a 6 

responding plaquette energies are given in column 4 of 
Table iTVl The energy for type 1 plaquettes is approxi- 
mately 50% larger (more negative) than for s = 2. In the 
extreme limit of very large and negative Vu the ground 
state has been determined previously in Ref|3. The state, 
referred to as the R state in Ref as well as in the remain- 
der of this manuscript, maximizes the fraction of type 1 
plaquettes, and is unique (up to lattice symmetries). 

The numerical investigation mentioned in the previ- 
ous subsection, shows that the R state is not the low- 
est energy state for the diagonal effective Hamiltonian 
at s = 3/2. Instead, we find a massively degenerate set 
of classical ground states. One example of these states 
has all the minority sites contained in a set of parallel 
Kagome layers of the pyrochlore lattice. Every Kagome 
plane will be have the same spin configuration shown in 
Fig. [21 This example, and many other states in this 
degenerate manifold all have a \fi x \/3 structure in the 
Kagome planes, and therefore we shall refer to a large 
subset of this manifold of states as the -s/3 X -s/3 states. 

The analysis of this degenerate manifold of states is 
somewhat involved. We therefore leave the details to 
Appendix |C1 and only mention a number of facts here. 
All the states we have found numerically have plaquette 
type fractions of xq = \,x\ — \-,%2 = %> x 3 = | an( i 
.i'4 = i. As a consequence, the energy per plaquette of 
these states is 

i^ x ^ = Uv 1+ 2V 2 + V^. (89) 
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FIG. 14: (Color online) 3:1 spin configuration of a single layer 
of tetrahedra in the \fS x ^/3 state. Only minority spin sites 
are marked by (red) solid circles. Flippable plaquettes (type 
1) are denoted by a (blue) circle drawn at their center. The 
same conventions as in Fig. 1121 are used here. 

In appendix|0 we show by explicit construction that the 
degeneracy is at least 

18 x 2t& + 4 x 3 L - 36, . (90) 

which grows exponentially with system size. We have 
not shown that the above states exhaust the possibilities 
with energy given by Eq. (|89|l . so the above formula is 
only a lower bound for the degeneracy. 



B. Effect of off-diagonal term 

In this subsection we add to the effective Hamiltonian 
the off-diagonal term where it is likely to be important 
(low values of s). For s = 1/2, the off-diagonal term is 
parametrically larger than the diagonal terms in the a <C 
1 limit. For s = 1, it is of the same order as the diagonal 
terms. However, our explicit calculations demonstrate 
that even in this case, the off-diagonal term is numerically 
more than four times larger than the largest diagonal 
plaquette energy. 

For s = 3/2, the off-diagonal term is negligible in the 
a <C 1 limit, but extrapolating the DPT results to the 
isotropic case a = 1 indicates that while it is not larger 
than the diagonal terms, it is likely not negligible either. 

To gauge the importance of the off-diagonal plaquette 
term, it is instructive to compare the diagonal energy of 
the various candidate ground states studied above. Their 
energies per plaquette are shown for small values of s in 
Table IVII We see that the energy differences amongst 
these competing states are rather small on the scale of 
the off-diagonal amplitude K. For instance, for s = 3/2, 
the energy difference between the "worst" of these three 
states (the R state) and the "best" (the x \/?> state) 
is only 0.0018 per plaquette, approximately four times 
smaller than the diagonal coupling K = 0.008. Thus we 
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TABLE VI: Diagonal energy per plaquette in various classical 
ground states. Energies are given in units of J z a 6 



can expect that adding the K term can introduce suffi- 
cient quantum fluctuations to alter the balance between 
these states, either destabilizing one in favor of the other, 
or perhaps stabilizing a superposition of these orders in 
some form. 

We cannot hope to establish the result of such sub- 
tle energetics here, particularly given the non-trivial na- 
ture of the effective QDM Hamiltonian including the off- 
diagonal term. However, we will discuss several natural 
candidate ground states from the perspective of order-by- 
disorder and the general theoretical framework of QDM- 
type models. 



Purely off-diagonal QDM - 



= 1/2 



Let us consider first the simple case of s — 1/2, for 
which the Hamiltonian is well-approximated by including 
the off-diagonal only. Clearly low-energy ground states 
of this Hamiltonian must have significant amplitude for 
type 1 plaquettes, as other plaquettes are annihilated by 
the off-diagonal term. We note that the trigonal state 
has no type 1 plaquettes. This implies that it is an ex- 
act zero energy eigenstate of the purely kinetic Hamilto- 
nian. Since it is straightforward to construct states with 
significantly negative energy per plaquette, the classical 
trigonal state is clearly an excited state in this case. It 
seems difficult to imagine any way that the ground state 
could be adiabatically connected to the trigonal state (or 
any other zero energy state with no type 1 plaquettes). 

Let us instead consider what sorts of states might nat- 
urally minimize the energy of the kinetic term. This sort 
of pure QDM problem has been considered in numerous 
places in the literature. Specifically for the QDM on the 
diamond lattice, the question has been discussed in Ref. 33 
(see references therein for a guide to QDMs). Roughly 
speaking, the energy is minimized by delocalizing the 
wavefunction as much as possible amongst different dimer 
configurations. However, the non-trivial connectivity in 
the constrained space of dimer coverings makes the na- 
ture of this derealization subtle. 

One possibility in such a 3d QDM is that the ground 
state is a U(l) spin liquid, in which the derealization is 
sufficiently complete as to prevent any symmetry break- 
ing (the meaning of the U (1) is discussed in-depth in e.g. 
Ref. 0) ■ Roughly speaking, the wavefunction has support 
for all possible dimer coverings, with equal amplitude for 
all topologically equivalent configurations. The existence 



and stability of such a state can be established in a QDM 
with a particular form of diagonal interaction, in the 
neighborhood of the so-called "Rokhsar-Kivelson" (RK) 
point. While this point (corresponding to Vi = K > 0, 
Va = V4 = 0) is not physically relevant to the pyrochlore 
antiferromagnets, it is possible that such a U{1) spin liq- 
uid state remains the ground state for the purely off- 
diagonal QDM. 

A second possibility is that the derealization is incom- 
plete, due to "order-by-disorder" physics. In particular, 
it may be favorable to delocalize only over a limited set 
of classical states, amongst which the connections are 
greater than those amongst generic classical configura- 
tions. In this case there is generally some symmetry- 
breaking induced by the selection of the states involved. 
Two sorts of such ordering have been proposed and ob- 
served in other similar QDM models. The first type of 
order-by-disorder state is one in which the set of classical 
states for which the ground state wavefunction has the 
largest amplitude are "centered" about a single classical 
state having the maximal number of type 1 plaquettes. 
Such a wavefunction may be "selected" by the kinetic 
energy, since under the action of the kinetic term of the 
QDM, this is the classical state is connected to the largest 
number of other classical states. In our problem, this 
classical state is just the R state mentioned above and 
discussed at length in RefslallOl A simple form for such 
a wavefunction is 



\R-Aip}) = cx p 



ip 



h.c. 



|R),(91) 



where |R) is the classical R state (with definite Sf = 
Sat), and 7p are variational parameters which can be 
used to optimize the quantum state |R, {jp}). 

The second type of order-by-disorder state is one in 
which there are a maximal number of independently res- 
onating plaquettes. This is based on the observation that 
the exact ground state for the kinetic term on a single 
plaquette is simply an equal amplitude superposition of 
the two type 1 states. However, neighboring plaquettes 
share sites, and therefore it is not possible to form a 
direct product of such resonances on all plaquettes. In- 
stead, the best one can naively do along these lines is to 
find the classical state with the largest number of type 
1 plaquettes which can be independently flipped, and on 
these type 1 plaquettes form an equal amplitude super- 
position of these two states. 

A state with the maximal number of independently 
flippable plaquettes can be described, starting from the 
V3 x V3 states introduced in Sec IV A3l The largest set 
of independently flippable plaquettes is a subset of all 
the type 1 plaquettes. An appropriate choice in a sin- 
gle plane is demonstrated in Fig. ^3 which includes half 
of the flippable plaquettes in the plane. It is interest- 
ing to point out that in each plane there are 2 possible 
choices of the plaquettes to be resonated (one half or 
the other), so out of each x state we can con- 
struct 2 L different choices of the plaquettes that will be 



resonating. The degeneracy of these state therefore is 
2 L x 3 L x 4 = 6 L x 4. Other states realizing this maximum 
number of independently resonating plaquettes may be 
possible, but we have not pursued this further. We refer 
to these states as "Resonating Plaquette States" (RPS). 
A precise wavefunction describing the RPSs we have de- 
rived from the V3 x y3 states is 

\RPS)= II -4(l + (|0^)(0 B |+h.c.))|*) , (92) 

PGG VZ 

where G denotes the set of non-overlapping resonating 
plaquettes, and \*S?) denotes one of the -s/3 x V3 states. 
There are 4 x 3 L choices for \^), and 2 L choices for G 
given | 'J). We note that the symmetry of the RPS is 
distinct and lower than that of the v^x v3 state - even in 
a single layer. Thus there is a precise distinction between 
these two states independent of the detailed form of their 
wavefunctions, for which the above explicit forms are of 
course only crude approximations. 

While potentially there might be some other state we 
have not anticipated, we think that most likely one of 
the three above states obtains in the purely kinetic QDM 
valid for s — 1/2. We will, however, refrain from mak- 
ing any definite statement as to which of these is the 
true ground state. One may imagine comparing the en- 
ergies of the wavefunctions in Eas. (|91l92|) to gauge the 
relative favorability of the R and RPS states. Unfortu- 
nately, even evaluating the variational energy of the |R) 
state in Eq.JMJ is rather challenging. Another difficulty 
is the considerable freedom in choosing the RPS wave- 
functions. Furthermore, a good variational wavefunction 
for the spin liquid is also needed for a more complete 
comparison. As always, there is much arbitrariness in 
defining each variational wavefunction, making the pre- 
dictive power of such an approach unclear. We believe 
this issue is more likely to reliably resolved in the future 
thorough numerically exact methods such as quantum 
Monte Carlo or exact diagonalization. 



2. s > 1/2 QDMs 

For s = 1 and s — 3/2, significant diagonal terms enter 
the QDM Hamiltonian. These act to alter the balance 
between the three candidate states discussed above, and 
also potentially to introduce the possibility of other states 
disfavored in the purely kinetic Hamiltonian. For both 
s = 1 and s — 3/2, the ground state of the classical 
diagonal term alone is actually a massively degenerate 
set of states discussed briefly in Section IV A 31 and in 
more detail in Appendix. IO From all the -s/3 x -s/3 states, 
we can construct RPSs. For s = 3/2, as we have shown, 
however, the R state is also quite low in diagonal energy 
and indeed only slightly worse than the x -\/3 states, 
as far as the diagonal term is concerned. Thus we expect 
that introducing the diagonal terms tends to favor both 
the RPS and the "renormalized" R state over the U(l) 



21 




FIG. 15: (Color online) Choice of non-overlapping flippable 
plaquettes to resonate in a plane of the \/3 x y/3 state. The 
chosen plaquettes are marked with (red) crosses in the middle. 



spin liquid. If their effects are strong enough, they could 
also stabilize the "non-resonating" x \/?> states (or 
any one of the other states with the same energy). We 
speculate that a spin liquid is unlikely to be realized in 
these cases, but that the RPS, R, and V3 x \/3 states (or 
more precisely all the states degenerate with the -\/3 x \/3 
states) remain very reasonable candidate ground states 
for these values of s in the isotropic limit. 



VI. DISCUSSION 

Since the development of the DPT and its analysis 
in this paper is rather involved, we begin in the first 
subsection by recapitulating the central points. In the 
second subsection we will then turn to a brief discussion 
of the implications on experiments and future directions 
of this work. 



A. Summary 

As a prototypical model of a magnetization plateau in 
a strongly frustrated quantum antiferromagnet, we con- 
sidered in this paper a nearest-neighbor spin-s model on 
the pyrochlore lattice at half the saturation magnetiza- 
tion. Such plateaus have been observed in the spinel 
materials HgCr 2 04 , and CdCr 2 04. We argued that a 
useful starting model is the easy-axis XXZ Heisenberg 
model in an external field, Eq. (0J . This model possesses 
all same symmetries as the isotropic Heisenberg model 
in an external field, and indeed we were able to extrapo- 
late our results to this limit. This model has the advan- 
tage that the transverse spin fluctuations can be treated 
systematically as a perturbation to the underlying Ising 
model. The resulting Ising model can be rewritten as a 
sum over the elementary tetrahedra of the pyrochlore lat- 
tice. In this Ising limit on the plateau, the spins on each 
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tetrahedron satisfy a 3:1 constraint, comprising a set of 3 
majority spins fully polarized parallel to the field, and 1 
minority spin antiparallel to the field. The half-polarized 
state has a macroscopic degeneracy corresponding to the 
number of possible positions for all the down pointing 
spins in the lattice. It is expected that the transverse 
spin fluctuations will play a role in selecting a ground 
state or set of ground states from the massively degener- 
ate 3:1 manifold. In this way, we are lead to a theoretical 
model involving a "constrained" degenerate perturbation 
theory in the 3:1 manifold. Our paper is devoted to a 
detailed analysis of such a theory and many parts are 
couched in sufficiently general language to be applicable 
to a broad class of systems. 

We began our discussion of the constrained easy-axis 
degenerate perturbation theory by deriving the general 
structure of the effective Hamiltonian that occurs at each 
order of dimensionless coupling a = J±/J z , Eq. I|15fl . We 
found that the effective Hamiltonian could be cast into 
a convenient form by performing a unitary transforma- 
tion that rotates all down pointing (minority) spins to 
up pointing spins and also by introducing a connectiv- 
ity matrix (whose elements are one for nearest neighbor 
spins and zero otherwise). The latter makes it possible 
to convert the sums over nearest-neighbor lattice sites to 
sums over the entire lattice (|23|) . These transformations 
cast the terms of the effective Hamiltonian coming from 
each order of perturbation theory into a form rather con- 
venient for analysis. The resulting terms are expressed 
explicitly in terms of the Ising variables on the lattice 
sites, the spin s, and the connectivity matrix. These 
terms were studied order-by-order in perturbation the- 
ory. We found that diagrams representing these terms 
naturally fell into two categories: contractible and non- 
contractible. Contractible diagrams are those whose de- 
pendence on some of the Ising variables is eliminated by 
summing with respect to their site index over all lattice 
points. Thus, a function of N Ising variables can be re- 
duced to a function of less than N Ising variables after 
this "contraction" process. The allowed contractions de- 
pend on the lattice geometry, the 3:1 constraint, and the 
Ising nature of the spin variables. Diagrams for which it 
is not possible to perform a contraction (equivalently, a 
reduction in the number of relevant Ising variables) we 
termed non-contractible. 

The central result of the analysis of contractible and 
non-contractible diagrams is that all terms in the con- 
strained degenerate perturbation theory up to and in- 
cluding 5th order are constant within the 3:1 manifold. 
Individual terms are shown to be constant by first con- 
tracting the diagrams as much as possible and then not- 
ing that the value of the diagram is unchanged under 
permutation among site indices associated with the Ising 
variables. The latter statement implies that the value 
of the diagram is independent of spin configurations al- 
lowed in the 3:1 manifold and hence a constant. In a 
similar manner, most terms at 6th order are also shown 
to be constant. However, we also observe that, at 6th 



order, there appears a "single large loop" diagram which 
cannot be contracted, and also defies the permutation 
arguments mentioned above. In fact, this loop diagram 
brings about non-constant contributions to the effective 
Hamiltonian in the 3:1 manifold. Therefore, this is the 
lowest order term which lifts the degeneracy of the 3:1 
manifold (at least for s > 1). The non-constant 6th or- 
der term includes effective interactions among spins on 
each hexagonal plaquette of the pyrochlore lattice. De- 
pending on the arrangement of minority sites, there are 
five distinct kinds of plaquettes that may appear and we 
label them 0,1. ..4 (See TablcQJ. Using the results of our 
degenerate perturbation theory, we evaluate the energy 
of each of these plaquettes as a function of a and s, and 
correct a mistake in Ref. The 3:1 condition constrains 
the allowed ratios of the various plaquettes in the lattice 
and allows us to express the total energy of the system 
(up to an overall constant), Eq. H6Uf) . in terms of only 3 
energies (1^21 . 

As a check on the results immediately above and as a 
further test of the robustness of those results, we also per- 
formed a large-s expansion in the easy-axis limit. As with 
the fully quantum theory, we expanded the harmonic spin 
wave energy in powers of a up to the 6th order, applied 
the diagrammatic analysis above involving contractible 
and non-contractible diagrams, and studied the resulting 
energy of the non-constant 6th order terms. The result 
(|80|l agrees exactly with the O(s) term obtained from 
the quantum degenerate perturbation theory, <|64|l . This 
satisfying consistency tells us that the large-s limit and 
small-a limit commute, and thus our analysis is likely 
well controlled. 

In the final section of the paper we used the results of 
the degenerate perturbation theory to determine the low 
energy states on the plateau as a function of s. Our result 
that the first non-constant diagonal term in perturbation 
theory comes at 6th order is independent of the spin value 
s. However, terms that allow plaquettes (such as type 1) 
to resonate occur at order 6s, which can be either larger 
or smaller than 6 depending on s. In the strict easy-axis 
limit, therefore, for s > 3/2, the low energy states are 
therefore determined only by a diagonal effective Hamil- 
tonian, which can be analyzed classically. In the large 
but finite s limit, we are able to resolve the degeneracy 
of the "zero flux" manifold found in the large-s analysis 
(extended from that of Hizi and Henley^ to the XXZ 
model). We predict a "trigonal 7 " state (see Fig. IT^ to 
be the exact ground state in this easy-axis limit and for 
large s, and numerical analysis suggests this obtains for 
s > 2. For s = 3/2, the lowest energy configuration we 
have found in the Ising limit is a massively degenerate set 
of states (for example the v^x -\/3 states, see FiglTlll. For 
s < 1, and for s = 3/2 extrapolated to the isotropic limit, 
we find that the off-diagonal term in the effective Hamil- 
tonian becomes significant, and we suggest several likely 
candidates for the ground states in these cases. This in- 
cludes a possible U(l) spin liquid state, which would be 
quite remarkable if realized. 
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B. Implications and future directions 

First let us comment briefly upon the relevance to 
the spinel chromites. For HgCr2 04, it is known that 
the temperature at which the plateau forms (w 7°K) is 
comparable to the highest temperature at which mag- 
netic order is observed. The theoretical estimate of the 
magnitude of the couplings in the effective Hamiltonian 
due to quantum fluctuations for s — 3/2 is however 
small, e.g. V x « 0.02 J from Eq.^. Thus the tem- 
perature at which quantum fluctuations are expected to 
induce magnetic ordering would be very low. A crude 
estimate based on the measured Curie- Weiss tempera- 
ture in HgC^O^ would predict an ordering temperature 
< 0.2K. This strong quantitative disagreement with ex- 
periment indicates that a stronger classical mechanism 
- i.e. physics outside the Heisenberg model - must be 
behind the plateau formation. Indeed, a recent study of 
a simple model of spin-lattice coupling gives a reasonable 
explanation of the plateau and its order, predicting sta- 
bilization of the R stated It would be quite interesting 
to see whether quantum fluctuations might however play 
a role in the other chromite spinels, e.g. CdCi'204. 

We now move away from the experiments on HgCr 2 04, 
where the pure nearest-neighbor Heisenberg antiferro- 
magnet neglecting spin-lattice interactions is clearly in- 
adequate. Instead, we would like to address a basic ques- 
tion that may be in the mind of the reader. In the pure 
spin-s isotropic Heisenberg model (i.e. J± = J z = J), 
is there a plateau at half-magnetization? At s — 00, i.e. 
the strict classical limit, the answer is no, and indeed the 
magnetization is a simple linear function of field in this 
case. In principle this question can be addressed by the 
1/s expansion. However, to the order studied, the situ- 
ation remains unclear: the leading-order spin- wave spec- 
trum remains gapless even in a field. Higher-order cal- 
culations in 1/s are required to resolve this question via 
that approach. Within the XXZ model, for any amount 
of anisotropy (a < 1), a plateau is expected even in the 
classical limit, so by continuity it is likely to persist at 
smaller s. However, the extrapolation to a = 1 is not 
clear. In Appendix ^ we present some simple calcula- 
tions aimed at addressing the plateau width. In particu- 
lar, we show that the plateau narrows both from above 
and below upon perturbing away from the Ising limit, 
where it is maximal. The plateau edges are determined 
by the points at which the gap to excitations with non- 
zero S z vanishes. Unfortunately, unlike the calculation of 
the splitting within the plateau states (the main focus of 
this paper), the energy difference between the plateau 
ground state and excited states with higher/lower S z 
is non- vanishing already at quadratic/linear order in a. 
Hence, a high-order calculation of this gap becomes much 
more involved than those in the bulk of this paper, and 
an extrapolation to the isotropic limit is probably not 
reliable. The existence of a plateau in the isotropic limit 
is a subject worthy of study by other methods. 

Next we turn to future applications of the formalism 



developed here to other problems. From our exposition, 
it should be evident that our methods generalize rather 
straightforwardly to other models of quantum antiferro- 
magnets with Ising anisotropy, provided a few conditions 
hold. First, the lattices should be composed of site- 
sharing simplexes. A simplex is a collection of sites in 
which every pair of sites is connected by a bond; exam- 
ples include triangle, crossed square, and tetrahedron. 
Second, the ground states of the Ising part of the Hamil- 
tonian on a single simplex should all be permutations 
of one another. This allows Ising exchange, single-site 
anisotropics, biquadratic and other interactions. Third, 
the interactions should be the same on each bond, but 
could include quite arbitrary combinations of exchange, 
biquadratic couplings etc. There are quite a number of 
interesting models of frustrated magnetism which share 
these features. For instance, the XXZ models on the 
Kagome and checkerboard lattices can be studied this 
way at several values of the magnetization. The XXZ 
model on the pyrochlore lattice at zero field is also such 
a system. It will be interesting to explore the behavior 
of these models at various values of s using the methods 
of this paper. 

More generally, the methods of this paper are possible 
because of a key simplification: in a strong magnetic field, 
the symmetry of the spin Hamiltonian is U(l) rather than 
SU(2). Many more theoretical methods are available to 
treat systems with abelian conserved charges than for 
5J7(2)-invariant spin models. Furthermore, in the inter- 
esting search for spin-liquid states of quantum antiferro- 
magnets, much theoretical success has been achieved in 
recent years in realizing such states in t/(l)-symmetric 
models, while examples of 5?7(2)-invariant spin liquids, 
even in models, are much more limited. Therefore it 
seems likely that quantized magnetization plateaux may 
be an excellent hunting ground for such exotic states of 
matter, and moreover there is hope for theory and exper- 
iment to meet on this plain. 
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APPENDIX A: EDGES OF THE 
MAGNETIZATION PLATEAU 

In this appendix we analyze the gap to spin excita- 
tions above the 3:1 plateau. To lowest order in the trans- 
verse spin fluctuations, in the field range corresponding 
to the plateau in the Ising model, excitations carrying 
AS Z = ±1 are are separated by an energy gap of or- 
der J z from the manifold of half-polarized states. As 
the magnetic field is varied this gap decreases. At zero 
temperature, the magnetization plateau ends when the 
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energy gap to one of these spinful excitations vanishes. 
We shall consider the spin AS = ±1 "single particle" ex- 
citations above the 3:1 manifold only, and find a rough 
perturbative estimate for the limits of the half magneti- 
zation plateau region. Clearly, the plateau region cannot 
extend beyond the magnetic field values at which the gap 
to these excitations vanishes, and so can only be more 
narrow than the extent we shall find here. 

Consider AS = +1 excitations above the ground state. 
These excitations are more favorable in higher magnetic 
fields. The simplest way to insert such an excitation is 
by raising one minority spin by one unit, i.e. obtained 
by acting with Sj on a minority site to change Sj = —s 
to Sj = 1 — s. This creates two "defective" tetrahedra 
(shared by this site), which no longer have the optimal 
3:1 spin configuration. We will call the space of such 
states "manifold A". For s > 1/2, manifold A comprises 
all the lowest energy AS = +1 excitations in the Ising 
Hamiltonian. Other AS = +1 states which involve more 
spin raising/lowering operators, e.g. those created by 
Sj~Sj~S^ , have higher energy because they either create 
more defect tetrahedra or make the two defect tetrahedra 
more energetically costly. The exception is s — 1/2, for 
which manifold A is actually incomplete, and there are 
other AS = +1 excitations outside it with the same th 
order energy. We will henceforth assume s > 1/2 in the 
remainder of this appendix to avoid this complication. 

The th order energy difference from the ground state 
manifold is 



E Q A -El x = 2J z {Zs-h) 



(Al) 



The gap shrinks as h increases, and vanishes at a mag- 
netic field Ha = 3s. This corresponds to the high field 
edge of the half polarized plateau in the a = limit. The 
single site excitation is however highly degenerate. The 
excitation can reside on any one of minority sites (N 
is the number of pyrochlore sites). We expect breaking 
this massive degeneracy (i.e. having this magnon "flat 
band" acquire some dispersion) will lower the energy of 
this excitation (closing the gap at lower magnetic field). 
One therefore needs to take into account corrections from 
higher orders in DPT. 

To first order in a, only spin ASf = ±1 can be trans- 
ferred from one site to a nearest neighbor. In any 3:1 con- 
figuration, the nearest minority spins reside on next near- 
est neighboring sites in the lattice sense, which means 
these sites are a distance of two links apart. Therefore 
this process always leaves manifold A, and the first order 
term has no contribution. 

To 2nd order in a there are various viable processes. 
Apart from hopping spin 1 between neighboring sites, 
and then back, there is also a non-trivial process in which 
the defect site (the — s + 1 spin) swaps with one of the 
(next-) nearest neighbor minority spins. This is also a 
legitimate member of the A manifold. Processes not in- 
volving the defect directly are affected by the presence 
of the defect only when involving the nearest neighbors 
of the defect site. All other sites have the same contri- 



bution to the 2nd order correction as in the ground state 
manifold. 

We immediately conclude the form of the effective 
"magnon Hamiltonian" to 2nd order is 



i_/2 i_/2 

H A - H 3:l - 



Cl 



+ C 2 J2 { e i e 3 + h - C ) 



, (A2) 



where ((ij)) denotes next nearest neighbor sites on the 
pyrochlore lattice that are both minority sites, C1.2 > 
are coefficients depending on the physical couplings, and 
the operator ej creates a local excitation - it replaces a 
spin Sj = —s with a spin Sj=l — s on site j: 



a 



(A3) 



An arduous yet straightforward calculation for the A 
manifold results in 



H-A - ™3:1 = 

J z a 2 \ 3s (4s 2 - 3s + l) 



2s- 1 



4s - 1 



c 2 , 



e\e + h.c 



(A4) 



Note the singularity in Eq. I|A4(1 for s = 1/2, which re- 
flects the incompleteness of manifold A in this case. 

In a hopping Hamiltonian of the form of Eq. (|A4|) , the 
magnon eigenstates are delocalized Bloch states. The 
spectrum of these states in general depends in detail upon 
the structure of the lattice of minority sites. However, if 
we are interested in the minimum energy state only, there 
is an amusing simplification. Since the hopping above is 
everywhere negative, the lowest energy state is expected 
to be nodeless, Thus the minimum energy state is simply 
a constant amplitude (i.e. zero quasimomcntum) Bloch 
state. For any configuration of the minority spins |^o) 
in the ground state manifold, the corresponding state in 
the A manifold is 



I*) 



Direct calculation yields 

£ (e\e J+ h.c) |*) = 6|tt), 



(A5) 



(A6) 



where 6 comes about as the number of next nearest neigh- 
bors with spin S* = —s the excitation can hop over to. 
Because in any 3:1 configuration there are always pre- 
cisely 6 minority spins two links away, we are therefore 
able to obtain the lowest-energy AS Z = +1 magnon en- 
ergy irrespective of the particular arrangement of the mi- 
nority sitesl 

The minimum energy of i|A4fl is therefore 



Ei 



E 3:l 



-J z a z 



. 3s 



4s + 1) 



(2s- l)(4s- 1) 



(A7) 
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The combined gap to manifold A is now 



AE A = J z 2(3s -h)- J z a- 



3s i 



-4s+l 



(2s- l)(4s- 1) 



The gap vanishes at a magnetic field of 



(A8) 



E B = —J z sac, with c > a dimensionless constant (the 
largest eigenvalue of Wy). This gives 



hs = s(l + ac) > 



(A13) 



3s — a 2 



3s (8s 2 - 4s + 1) 
2(2s- l)(4s- 1) 



(A9) 



Note that for s ^> 1, the 0(a 2 ) correction is well- 
behaved, i.e. it scales in the same way as the zeroth 
order threshold field. 

Excitations with AS = — 1 (manifold B) can be real- 
ized by replacing a spin S*| = +s with a spin Sj = s — 1. 
The Oth order energy difference from the ground state 
manifold is calculated in an identical manner as in the 
AS = +1 case, and yields 



Thus the upper edge of the plateau decreases and the 
lower edge of the plateau increases as a is increased, nar- 
rowing the plateau with increasing quantum fluctuations. 
It is plausible that for sufficiently large a the plateau is 
obliterated entirely, the upper and lower edge meeting. 
Unfortunately, this is beyond the realm of this pertur- 
bative analysis. Even a naive extrapolation of the above 
lowest-order results is non-trivial, since the lower plateau 
edge depends upon the non-trivial constant c. It can in 
principle be computed for the various states obtained in 
Sec El but we do not do so here. 



h B ~ ^3:1 



2J z (h-s) 



(A10) 



At this order, the AS = ±1 excitation spectra are sym- 
metric about h — 2s. 

Oppositely from the behavior at the high field edge of 
the plateau, the gap shrinks as h decreases, and vanishes 
at a magnetic field hs = s. Once again, due to a huge 
degeneracy ( the defect can reside on any one of ^ ma- 
jority sites) we must resort to DPT in order to break the 
massive degeneracy of manifold B, expecting to lower the 
excitation energy. 

In contrast to the AS Z = +1 excitations above, it is 
already possible at first order to hop the spin S z = s — 1 
onto other sites and stay within the B manifold. This 
is because the majority sites are not isolated (4 of the 
6 neighboring sites of a majority site are also majority 
sites). Therefore, for any given 3:1 state, one can imme- 
diately obtain the effective Hamiltonian to 0(a): 



(All) 



Here Wij is the connectivity matrix of the lattice of ma- 
jority sites for the particular 3:1 state in consideration, 
and the operator creates a local excitation - it replaces 



a spin S z = s with spin Sj = s — 1 on site j: 

h) = \ S -l)j{s\j. 



(A12) 



Since a > 0, the hopping amplitudes in this case are 
positive, rather than negative as above. Thus, unfortu- 
nately, the lowest-energy eigenstate is a non-trivial Bloch 
state, whose energy depends upon the precise form of 
Wij, i.e. it is different for each of the 3:1 states. There- 
fore it is difficult to say anything specific about the 0(a) 
correction to the minimum AS" 2 = — 1 magnon energy. 
Because one can easily construct variational states with 
negative hopping energy (e.g. an antibonding state with 
support only on two sites), the 0(a) correction must 
be negative, and clearly from Eq. (|A1 l|l is of the form 



APPENDIX B: ALTERNATIVE CALCULATION 
OF THE DEGENERATE PERTURBATION 
THEORY 



In this appendix we present an alternative way of cal- 
culating the DPT 6th order diagonal term, providing a 
check on the calculation described in Section ITTTl 



The perturbation © can be viewed as a sum over spin 
transfer operators 



Si Sa 



h.c 



Mm") 



hi + h.c. 



(Bl) 

where hp transfers one quantum of spin angular momen- 
tum from site i to site j if these are on adjacent sites 
i,j. As a result, the n-th order virtual processes can be 
classified by choosing n links to act on with spin transfer 
operators, in a particular order. In order to return to 
the initial 3:1 configuration, the links must form closed 
loops of spin transfer. These include self retracting loops, 
defined as loops where after following some path in the 
lattice, we turn back and travel the same path in the re- 
verse order back to the origin. The chosen links can be 
represented by a graph on the lattice, by coloring these 
links. Apart from the graph, we must also specify the 
order by which the links operate. 

The virtual states in a DPT process will have a number 
of sites in a different spin state relative to the initial 3:1 
state. We denote Sj — Ujis—rrij), where the variables rrij 
are non-negative, and take on any integer value between 
and 2s. Let us denote the set of modified sites by 
M . Thus rrij ^ only if j 6 M. Using the fact that 

J2j = \ J2j 1 and a i a i = 0' bot]l due to tn e 3:1 

constraint, we find the inverse resolvent can be rewritten 
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FIG. 17: (Color online) Third order process. 



FIG. 16: (Color online) Second order process. 



as a sum involving only the modified sites 
Wo - E Q = J z ^2 SfSj - E' 

= Jz <JiO-j{s - mi)(s - TTLj ) - E' Q 

m 

— J z \ ^2 11 : n •">:"> j + 2s TOj j - Eg , 

\<ij>eM j'6Af / 

(B2) 

where we have absorbed all constant energies into E'q. 
With all rrij = 0, this resulting energy should vanish 
(energy difference to a state in the 3:1 manifold) and 
therefore conclude that the energy Eq = in the final 
expression above. Finally, our virtual state energy is 



Wo — E — J z 




2s ^2 m j 



(B3) 



At 1st order, the perturbation © takes any initial 
state out of the 3:1 manifold, and therefore gives no con- 
tribution (equivalently no loop can be closed with only 
one link). 

In 2nd order, we must act on the same link twice to 
undo the spin transfer, and return to the 3:1 manifold. 
Therefore all non-vanishing processes are confined to one 
tetrahedron. The same process can act on any one of 
the links of this tetrahedron, with the same resolvents 
for a given order of the link operators. Therefore, by 
summing up these diagrams, we will get a function of 
the spins on this tetrahedron which is invariant under 
any permutation of the 4 sites (the sum is represented 
in a diagrammatic way drawn in Fig. I16[) . Therefore the 
3:1 configuration on any single tetrahedron has the same 
energy at this order in DPT, no matter which corner is 
chosen to be the minority site. For this reason, the 2nd 
order contribution gives a constant shift in energy to all 
the 3:1 states. 

In 3rd order, the only way one can return to the 3:1 
manifold is by forming a single closed loop of spin trans- 
fer, around a triangular side of a single tetrahedron. As 
in the 2nd order processes, we can sum over all such pro- 
cesses occurring on the single tetrahedron. One again we 



end up with a function that is invariant under any per- 
mutation of the 4 sites of the single tetrahedron (the sum 
is represented in a diagrammatic way drawn in Fig. I17|l . 
Once again, we cannot distinguish energetically between 
the different 3:1 states defined on this single tetrahedron. 
As a result, the 3rd order term must produce a constant 
shift in energy. 

As exemplified in the 2nd and 3rd order terms, all DPT 
processes (of any order) that are confined to one tetrahe- 
dron (meaning all the spin transfer operators act on links 
in the same tetrahedron) all add up to constant shifts in 
energy. In every such case the argument is the same as 
above. Starting with a given process defined on a single 
tetrahedron, sum over all the processes of the same struc- 
ture on the same tetrahedron. Then one always arrives 
at a function that is invariant under any permutation of 
the 4 corners. Within the 3:1 manifold, these processes 
cannot favor one configuration over the other, and must 
produce a mere shift in the overall energy. More explic- 
itly, every process results in a function of the four Ising 
variables on a tetrahedron, as defined in Section [H] The 
sum of all diagrams with same structure then takes on 
the generic form 

/(fi> °~2, o"3i + ( a U permutations of 1...4) 
= a + ai (at + cr 2 + 03 + 04) 

+ a 2 (0J.02 + 0103 + (71(74 + (72(73 + O204 + O3C4) 
+ a 3 (71(72(73 (7 4 (<7i + (72 + <7 3 + (7 4 ) + a 4 (7l(72(73(74 • 

(B4) 

The 3:1 constraint gives ci 02 03 04 = —1 as well as 
(01 + 2 + 03 + 04) = 2. The generic 4 spin function 
then reduces to 

4 

/(01, 02, 03,04) = const + a 2 ^2 °i (J i 

«j=i 

2 



= const -\ — - j ^2 u 



const 



(B5) 



Given this result, we can always ignore virtual processes 
confined to one tetrahedron for the remainder of our dis- 
cussion. 

Similar arguments can be applied to processes confined 
to 2 adjacent tetrahedra. Summing over all possible pro- 
cesses confined to pairs of adjacent tetrahedra, one al- 
ways ends up with a constant energy shift in the 3:1 man- 
ifold. A pair of adjacent tetrahedra can only be in one 
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FIG. 18: (Color online) 4th order process. FIG - 19: (Color online) Sum of 4th order processes. 



of 3 configurations in the 3:1 manifold (see Fig. |3J). The 
adjacent tetrahedron pairs can be enumerated by the site 
connecting them. Therefore the number of distinct pairs 
is equal to the number of pyrochlore sites. In the 3:1 
manifold, the number of configurations where the shared 
site is a minority site (see Fig. |3(a)| is fixed to be \ of the 

the shared site is a 
with a 



of the sites, 



sites. In the remaining | 
majority site. The two configurations in Fig 
majority site connecting the pair of tetrahedra are distin- 
guishable physically because of the different directions of 
the links, but in terms of the diagrams defining DPT pro- 
cesses Fig. 3(b) and Fig. 3(c) are indistinguishable (they 
are identical in the graph sense - the link structure is the 
same). Therefore, the summation over all possible pro- 
cesses confined to these two adjacent tetrahedra always 
results in a function of the Ising variable defined on the 
shared site 



f(oi) = a + aio l . 



(B6) 



For each pair of adjacent tetrahedra, we have the same 
function JB6|) of the shared site Ising variable. Summing 
over all pairs of adjacent tetrahedra, is equivalent to sum- 
ming over the sites of the pyrochlore lattice. Therefore, 
the sum of all the DPT processes of this sort results in 

J2f(<T i ) = Na + a 1 Y, a i = Na o + Y ai > ( B7 ) 

i i 

where N denotes the total number of pyrochlore lattice 
points. Once again, we are led to the conclusion that all 
such processes can only give a constant shift in energy, 
and we shall ignore all such instances in the remainder 
of our discussion. 

At 4th order we have a number of possible processes. 
The only spin transfer processes comprising a single loop 
will be confined to 2 adjacent tetrahedra. According 
to the arguments above, this single loop process will 
give a constant shift in energy, within the 3:1 manifold. 
Apart from single closed loop, we can only also have pro- 
cesses with 2 separate closed loops, comprising two non- 
overlapping links which get acted on twice. Each chosen 
link occupies its own tetrahedron (if the two links are on 
the same tetrahedron, one must act on a bond with two 
majority spins, and automatically produces zero). The 
case in which these 2 tetrahedra are adjacent, will result 
in a constant, by our arguments above. Otherwise, there 
are only two possibility. In one case these two tetrahedra 



share a neighboring tetrahedron. In the 2nd case, they 
are separated by more than one tetrahedra. This classifi- 
cation is necessary because, in the former case, different 
resolvents may show up through the interaction term of 
the virtual state energy l|B3|) . If the two links are sepa- 
rated by more than one tetrahedron, this does not occur. 

In the case where the two tetrahedra are sufficiently 
well separated, with no mutual interaction in the virtual 
state energy, we can once again invoke the trick of sum- 
ming over all such processes on each single tetrahedron 
containing a link, and the DPT term cannot distinguish 
between the four 3:1 configurations on each tetrahedron. 
Once again we end up with a constant shift in the overall 
energy for all the 3:1 configurations. 

Now we consider the processes where the two tetrahe- 
dra share a neighboring tetrahedron. One such diagram 
is pictured in Fig. 1181 The structure spans a chain of 
three tetrahedra. In the present case, when summing 
over diagrams one must be cautious to sum over dia- 
grams with the same resolvents. Summing over equiva- 
lent diagrams on the chain of three tetrahedra, such as 
those found by permuting site 2 with 3 and 4, yields a 
function that treats the 3 sites at the outer side of each 
edge tetrahedron on equal footing. We graphically indi- 
cate the sum of diagrams, by coloring all the links that 
appeared in one of the diagrams we summed over (see 
Fig.HU). 

Formally, the sum of processes gives rise to a func- 
tion of the 8 Ising variables on the two edge tetrahedra 
/(oi, tj%). Because of the above arguments, this func- 
tion has to be invariant under permutations of the 3 outer 
edge spins. Considering all possible expressions we can 
form out of the 3 edge spins in a "symmetric" way 

/(cti , . . . , o 8 ) = a + a x (cr 2 + o 3 + o 4 ) 

+ a2(0"2<73 + 0204 + 0304) + G 3 O2O 3 '4, 

(B8) 

where the dependences of 00,1.2,3 on 0-1,5,6,7,8 are im- 
plicit. Using the 3:1 constraint on this tetrahedron 
(01 + 02 + 03 + 04) = 2, as well as the identities 
it implies CT1O2O3O4 = —1 (or 02O3O4 = —a{) and 
(02O3 + CT2CT4 + CT3CT4) = 1 — 2oi, we find that /(01, . . . , o 8 ) 
can be rewritten as a function of 01,5,6,7,8 alone; 

/(01, . . . ,o 8 ) = a + ai(2 - 01) + a 2 (l - 2o x ) - a 3 oi 
= &o + Si ci- 

(B9) 
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FIG. 20: (Color online) Summing over equivalent edges of a 
diagram to contract the dependence on sites in the cluster. 




FIG. 21: (Color online) Summing over equivalent edges of a 
diagram to contract the dependence on sites in the cluster. 



Repeating the same simplification for the three Ising vari- 
ables of the other edge, (76,7,8, one finds that the sum of 
diagrams above can only produce a function depending 
on the 2 spins connecting the edge tetrahedra to the third 
tetrahedron in the middle of the chain (sites 1 and 5 in 
Fig. I19fl . This function is therefore a function of 2 neigh- 
boring spins on a single tetrahedron /(ci, 05). Each three 
tetrahedron chain process can be represented by the py- 
rochlore link between the two sites shared between the 
tetrahedra ( sites 1 and 5 in Fig. I19|) . Summing over all 
the processes represented by the 6 links in the central 
tetrahedron which contains sites 1 and 5, one ends up 
with a function which is symmetric under an permuta- 
tion of the 4 sites on this tetrahedron. Repeating the 
steps outlined in Eqs. I|B4|) and 1B5(I . we find the result- 
ing function is constant. 

The particular case analyzed above is only an example 
of a vastly more general case, which we now describe. 
Given a set of diagrams that reside on chains of tetrahe- 
dra, with the two edge tetrahedra not sharing a neighbor- 
ing tetrahedron between them in any one of the diagrams 
(except through the chain itself), one can sum over sites 
at the edge that we are free to permute (see examples in 
Fig. 1201 and Fig. I21|) . The summation in both examples 
will result in a function of the spins at site 1. This pro- 
cedure of contracting the diagram can be continued from 
both edges, and we will always end up with a function of 
2 sites residing on one tetrahedron. We then sum over all 
the different choices of the last 2 remaining sites on this 
tetrahedron. We then end up with a symmetric function 
of the 4 sites on a single tetrahedron, which is always a 
constant, from Eqs. (|B4|I and (|B5|) . Therefore, we con- 
clude that all such diagrams can be summed over to give 
constant shifts in energy. We shall refer to these sets of 
diagrams as "retractable chains" . 

Bearing the above arguments in mind, we now turn 
to 5th order in DPT. Any diagram which consists of 
only one closed loop of spin transfer, is confined to two 
adjacent tetrahedra. As elaborated above, the different 
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FIG. 22: (Color online) 5th order process. 



3:1 configurations on the pair of adjacent tetrahedra are 
indistinguishable in the graph sense. These diagrams 
can therefore contribute only a constant shift in energy 
Other than these diagrams, the only possibility is to have 
two separate closed loops of spin transfer. One must in- 
volve a single link being acted on twice, and the other 
must be a 3-link loop. Both loops reside on two different 
tetrahedra. If the two tetrahedra are adjacent, or share a 
neighboring tetrahedron (as in in Fig. 1221) . the tetrahedra 
are part of a retractable chain. Our general analysis for 
retractable chains then applies here as well. Therefore, 
these diagrams must also sum up to a constant in the 3:1 
manifold. Otherwise the two tetrahedra are sufficiently 
well separated, and we can sum over the various diagrams 
on each tetrahedron separately. The summation gives an 
expression in which each of the 4 tetrahedron corners are 
on equal footing, and cannot distinguish between differ- 
ent 3:1 states, resulting in a constant energy shift. 

In agreement with the results of Section^ DPT terms 
of order less than 6 will not split the energy of the 3:1 
states. Finally, we turn to 6th order in DPT 

H 6 = -V (HiTZ) 5 H!V (BIO) 

where 1Z = 9,~ V 1 is the resolvent. At this order we have a 
new class of single closed loops - loops around hexagonal 
plaquettes of the pyrochlore lattice (see Fig. [2*3^1 , which 
will split the energy of different 3:1 configurations. These 
processes can be enumerated by the hexagonal plaque- 
ttes they act on. So we can write the contribution from 
these processes J2v9( ai • • • a& ) where the spins around 
the hexagonal plaquette are denoted 1 ... 6. 

All other single closed loop diagrams are confined to 
1, 2 or 3 adjacent tetrahedra, the first two we already 
know will sum up to constant shifts in energy. The self- 
retracting loop residing on a chain of 3 tetrahedra is de- 
picted in Fig. [2] and is also summed over to produce a 
constant, since this is a retractable chain. 

Next we consider the diagrams which comprise 2 closed 
loops. There must be two of length 3, or one of length 
4 and one of length 2. In the case of 2 loops of length 
3, each loop must reside on a single tetrahedron. If the 
two tetrahedra are adjacent, or even share a neighboring 
tetrahedron (as in Fig. l25fl . they form a retractable chain, 
and result in a constant. Otherwise the tetrahedra are 
sufficiently well separated so that we can sum over dia- 
grams on the two tetrahedra separately, and end up with 
a constant. 
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FIG. 23: (Color online) Hexagonal plaquette loop process. 
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FIG. 24: (Color online) Single self retracting loop on a chain 
of three tetrahedra. 



We now turn our attention to the case of one loop of 
length 4 and one of length 2. The loop of length 2 must 
be one link appearing twice, and is therefore confined to 
a single tetrahedron. The loop of length 4 must be con- 
fined to 2 adjacent tetrahedra (or even just 1). If the two 
loops are sufficiently well separated (by more than one 
tetrahedra), we can sum over all 2-loop diagrams in the 
tetrahedron where the loop of length 2 resides, result- 
ing in a constant. We can then sum over diagrams on 
the tetrahedron cluster containing the 4-loop to produce 
another constant. If the two clusters are adjacent, or 
separated by just one tetrahedron, we have a retractable 
chain, which we have already shown must result in a con- 
stant. 

The last set of 6th order diagrams are those of 3 closed 
loops, each comprising a single link, acted on twice. Each 
loop is therefore confined to a single tetrahedron. Us- 
ing the same arguments as before, if all three tetrahedra 
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FIG. 25: (Color online) 6th order diagram with 2 loops of 
length 3. 




FIG. 26: (Color online) 3 loop diagram residing on a closed 
chain of tetrahedra, enclosing a hexagonal plaquette. This is 
the "f" process. 




FIG. 27: (Color online) Diagram contributing to the function 



are separated by more than one tetrahedra from one an- 
other, we can sum over all such diagrams contained in 
the same tetrahedra. In particular, sum over diagrams 
related by permuting over sites on each single tetrahedron 
separately results in a constant. Even if only one tetra- 
hedron is well separated from the other two, we can first 
sum over diagrams on the isolated tetrahedron to produce 
a constant, and then deal with the other two tetrahedra. 
The remaining two will reside on a retractable chain of 
tetrahedra, again resulting in a constant. 

The only cases which we must deal with more carefully 
are those diagrams where each pair of tetrahedra either 
shares a neighboring tetrahedron, or the two are adja- 
cent. When at least one pair of loops resides on adjacent 
tetrahedra, the cluster of tetrahedra will always be a re- 
tractable chain, with either 4 or 3 tetrahedra. We are left 
only with diagrams where all three tetrahedra containing 
the loops are not adjacent, but rather have shared neigh- 
boring tetrahedra. There are 3 tetrahedron clusters pos- 
sible, where this occurs. Fig.l26lshows a diagram residing 
on a closed chain of 6 tetrahedra, enclosing a hexagonal 
plaquette. This is not a retractable chain. In Fig. 1271 the 
three tetrahedra containing the loops all share one sin- 
gle neighboring tetrahedron. Finally, in Fig. Elthe three 
loops reside on a linear chain of 5 tetrahedra. However, 
it is important to note at this point that the tetrahedra 
at the two edges of this chain (1 and 7) can be identified 
to give the cluster shown in Fig. [5U 

The processes on the cluster shown in Fig. cannot 
be contracted, and will split the energy of the different 
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FIG. 28: (Color online) Diagram contributing to the function 
"r". The loops are embedded in the tetrahedra marked 2,4, 
and 6. 



3:1 states. We denote the sum of processes with the 3 
loops on the same tetrahedra as in Fig. |2HI • • • 
There is also a contribution from all the processes where 
the 3 loops reside on the complementary set of tetrahe- 
dra, which must result in /((Tq, cr% ... £75) (the same func- 
tion with the spins cyclically permuted once). Cyclically 
permuting the spins again, takes us back to the set of 
processes we summed over in ]{o\ . . . o~q), and therefore 
/ must be invariant under two permutations. This is a 
good check for the correctness of the result we find. The 
sum f(ai ... erg) + f{&(n o~i . . . 05) accounts for all the di- 
agrams on a given cluster enclosing a plaquette, and the 
total contribution to the effective Hamiltonian can be 
written 



^2 ... ere) + /(c 6 , CTi . . . fJ 5 ) 



(Bll) 



The processes on the cluster shown in Fig. |^| can be 
contracted from the "loose ends" of the diagram, simi- 
lar to the contraction we have implemented for the re- 
tractable chains. For the particular example in Fig. 1271 
we sum over the corresponding diagrams with 2 permuted 
with 3 and 4. The resulting sum can only be a function of 
site 1, 5, 7, 9, and 12. We then sum over the correspond- 
ing diagrams with 7 permuted with 6 and 8, and then 
sum over 12 permuted with 10 and 11. After these sum- 
mations, the resulting function can only be a function of 
the sites 1,5 and 9. Finally, we sum over all choices of 
three such sites on the central tetrahedron, which must 
result in a constant (a symmetric function of the 4 sites 
of a single tetrahedron). 

Next, we turn to the processes on the chain in Fig. 1281 
At first glance, it would seem that these are retractable 
diagrams, that will amount to a constant. Such a summa- 
tion would include 81 different diagrams (we retract the 
tetrahedra 2,3,5, and 6, each time summing over three 
permutations of sites). However, two of these diagrams 
are of the form of Fig. [23 in which case the tetrahedra 
1 and 7 at the edges of the chain are identified. As a 
result, these two diagrams may have different resolvents 
from the other 79 diagrams in this sum, since the sites on 
tetrahedra 1 and 6 may now interact. These 2 terms have 
in fact already been taken into account in f{<T\ . . .ere), 
and therefore should not be added again to the effective 
Hamiltonian. The summation that makes contractible 
diagrams give constant energy shifts requires that the 
resolvents be identical for the entire set of diagrams. 

Rather than explicitly calculating the remaining 79 di- 
agrams, we employ the following trick. In order to pro- 



duce a constant, we add to the 79 diagrams two terms 
having the same resolvent, but replacing the spin vari- 
ables <7j with those at the positions 1 through 6 in Fig. 1261 
All these 81 terms have the same resolvent as that of the 
process in Fig. 1281 and therefore will result in a constant. 
The two additional terms do not represent actual DPT 
processes, and therefore must be subtracted from the ef- 
fective Hamiltonian. The terms we need to subtract are 
enumerated by the choice of central tetrahedron in the 
chain in Fig. 1281 (tetrahedron number 4), which we can 
identify with the tetrahedron containing sites 3 and 4 in 
Fig. [23 Given a choice of the central tetrahedron, each 
term of this sort is a function of the 6 spins r(u\ . . . <jq) 
as defined in Fig. [213 Given a hexagonal plaquette there 
are 6 choices of the central tetrahedron, and correspond- 
ingly 6 terms and we must take into account all of these. 
Choosing a different central tetrahedron is equivalent to 
cyclically permuting the spins 1 ... 6, and so the total 
contribution of all these process is 



v \ 



cti . . . er 6 ) + r(a 6 ... 0-5) 
+ r(<7 5 . . .04) + r(<7 4 . . .03) 
+ r(er 3 . . . a 2 ) + r(a 2 . . .01) 



(B12) 



The 6th order term in the effective Hamiltonian finally 
reads 

He = ^2 9(0-1 ■■■o-q) 
v 

+ ^2 (f( ai ■■■°~e)+ /(f 6 , CTi . . . cr 5 )) 

V 

- ( r ( CTl ■■■ <J e)+ r ( <T 6 ■■■0-5) 1 L> ' ' ' 



+ r(<7 5 . . . cr 4 ) + r(cr 4 . . . 03) 
+ r(a 3 . . . £7 2 ) + K°2 • • • 01) 



Calculating these functions explicitly, we find the exact 
same results as in Section [Hi 



APPENDIX C: SPIN s 



DIAGONAL TERM 



GROUND STATE DEGENERACY. 

In this appendix we analyze the lowest energy states 
found for the diagonal term in the effective Hamiltonian 
for spin s = § discussed briefly in Section IV A 31 The 
lowest energy states turn out to be massively degenerate. 
We have found what is at the very least a subset of this 
manifold, which already exhibits a degeneracy that grows 
with the system size, diverging in the thermodynamic 
limit. 

First we consider all the possible states we can con- 
struct with all the Kagome layers in the pyrochlore taking 



31 



on the V3 x V3 configuration in Fig. ^] Given a plane 
in this V3 x \/3 configuration, there are 3 different ways 
to place the following layer (also in the -\/3 x \/3 config- 
uration) above it, as described in Fig. I2H1 This freedom 
in the way the planar configurations are stacked is the 
source of a massive degeneracy - there are 3 L possible 
stacking choices, where L is the linear dimension of the 
system. Also, there are 4 plaquette directions in which 
to choose the direction of the stacking, resulting in an 
overall number of such states 4 x 3 L . For now, we shall 
work with this subset of the entire degenerate manifold, 
since these states are rather easy to handle. We shall 
refer to his set of states as the \/3 x -\/3 states. 

Now we turn to calculate the energy of these states, 
proving these are degenerate states. Each plaquette in 
the \/3 x a/3 planes shares links with three plaquettes 
above it. These 4 plaquettes enclose an up-headed cell, 
as defined in Section IV Al For any one of the 3 V3 x \/3 
spin configurations of the next Kagome plane, we always 
get the same up-headed cell types above the type 1 and 
type plaquettes in the V3 x V3 plane. The clusters 
enclosing these up-headed cells above the type 1 and type 



plaquettes are depicted in Fig. 29(a) and Fig. 29(b)| 
respectively. Since we are only considering up-headed 
cells, every one of the plaquettes appears in only one 
cell. 



Inspection of Fig. 29(a) shows that it includes one pla- 
quette of type 1, one of type 2, one of type 3, and one of 
type 4 (This i s just t he type 1 cell in Table |ffij . For the 
cluster in Fig. 29(b)| we find it includes two type pla- 
quettes, and two type 2 (type 8 cell in Tabic ITTT|) . Of the 
plaquettes in the planes, which comprise i of the plaque- 
ttes in the lattice, | are in a type 1 configuration, and ^ 
are in a type configuration. Since each planar plaquette 
determines the configuration of the unique up- headed cell 
it is a part of, we find therefore, that | of the up-headed 
cells are type 1 cells (yi = |) and the remaining i are 
type 8 cells (y$ = §)■ The corresponding plaquette type 
fractions are xq = |, X\ = \,xi = |, X3 = | and x± = ^. 
Finally, the energy per plaquette of the -\/3 x states 



If 



2V 2 + V A ) . 



(CI) 



We expect the fractions of plaquettes in all the states 
in this degenerate manifold to be the same as in this 
subset, since otherwise, it is rather unlikely (though not 
impossible) that a different combination of fractions will 
yield the same energy. 

Having analyzed this set of states, we turn to an ad- 
ditional set of states, with the same energy. Consider 
a particular subset of -s/3 x V3 states, with every two 
Kagome planes stacked in the same manner. This subset 
of the \/3 x \/3 states has a total of 4 x 3 2 = 36 states 
- factor 4 for choosing the direction of the layering, a 
factor 3 for the choice of how to position the planar con- 
figuration on one plane, and another factor of 3 from the 



freedom to choose how to stack the next planar layer. 
We shall refer to these as the uniformly stacked states. 

Starting from any one of these uniformly stacked states 
we note that the type 1 cell can swap its type 1 and 
type 3 plaquettes by changing the position of only one 
minority site. In Fig. 29(a)| this can be accomplished by 
moving the minority site at site 4 to site 5. In addition, 
site 5 is part of a type 8 cell, and it is denoted also in 
Fig. |29(b)[ where it is evident that in order to maintain 
the 3:1 constraint, wc must also shift the minority site 
at site 1 in Fig. 29(b) to site 6. Because the stacking of 
the next layer is exactly the same, the same shifting of 
minority sites must occur along the entire straight line 
passing through site 4 in the direction from site 4 to site 
5. Closing this chain at infinity makes this chain into an 
infinite length loop of alternating minority and majority 
sites, that are flipped, and therefore this maintains the 
3:1 constraint. One can convince one self from Fig. 29(b) 
that the type 8 cell will remain a type 8 cell under these 
minority site shifts. Since all the cells remain in the same 
configuration type, the plaquette type fractions remain 
the same as in the x states, have the same energy, 
and are therefore degenerate. 

We now turn to calculate the degeneracy of this new 
set of states. Each type 8 cell has exactly one site it can 
shift in this manner, but since the chains are not shared 
between different type plaquettes in the plane, one can 
convince oneself that these straight line chains of sites 
can be flipped independently. Therefore, starting from 
a particular uniformly stacked state, since there are 
plaquettes in every plane (where N is the number of py- 
rochlore sites, and the number of hexagonal plaquettes), 
there are j^j- chains we can independently flip between 

two configurations, resulting in a degeneracy of 36 x 2 Tar . 
However, an additional subtlety must be addressed. 

Every V3 x -\/3 state has all the up pointing tetra- 
hedra in one of three 3:1 configurations. Flipping any 
chain in a uniformly stacked state will introduce some 
number of up pointing tetrahedra in the 4th possible 3:1 
configuration of a single tetrahedron. Starting from one 
uniformly stacked state, flipping all possible chains, we 
will change all the tetrahedra of one of the three 3:1 con- 
figurations into the 4th tetrahedron configuration type, 
which was absent in the initial uniformly stacked state. 
Once again we will find ourselves with all tetrahedra in 
only three possible different 3:1 configurations. This sug- 
gests, that perhaps this final state we have reached is also 
a uniformly stacked state. From inspection we find this 
indeed is the case, so that from one uniformly stacked 
state we can reach one other such state by flipping all 
possible chains in the manner described above. We must 

N 

therefore correct the degeneracy to 18 x 2t5t to account 
for this double counting. 

In total we have found the degeneracy of these two sets 
of states to be 



18 x 2T5r +4X3 1 



36 , 



(C2) 



where we have subtracted 36 since this is the number of 
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(a) Stacking choices (b) Stacking choices 

above a type 1 above a type 

plaquette. plaquette. 



FIG. 29: (Color online) Stacking of two planar \/3x v3 states. 
The solid lines represent the tetrahedra of the lower plane. 
The dashed lines represent the tetrahedra of the upper layer. 
The same convention regarding up pointing and down point- 
ing tetrahedra applies here as in Fig. [H for both the solid 
and dashed line triangles. In both figures, the dashed tetra- 
hedron with the corners 1, 2,3 must have one of these three 
corners a minority site, and any one of these three can be cho- 
sen. Once the minority site has been chosen from 1, 2, 3, the 
\/3 x %/3 state of the upper layer is then uniquely determined. 
The choice of site 1 is explicitly shown in both cases. 



states that appear in both sets of states we have analyzed 



(these are simply the uniformly stacked states). We shall 
refer to the combined set of states as the 1 — 8 manifold 
of states, since it involves only type 1 and type 8 cells. 
The degeneracy we have calculated matches precisely the 
degeneracies in Table IV1 for all the clusters where we have 
found the lowest energy state for s = |, namely the 3 x 
3x1, 6x3x1, 3x3x2 and 3x3x3 clusters. Our 
analysis clearly shows that this set of states is massively 
degenerate, and despite exhausting all the states we have 
found numerically with this energy, we cannot be certain 
that these exhaust all possible states with this energy. 

In a previous publication^ analyzing the same 3:1 de- 
generate manifold of states, the authors ascertained the 
maximum fraction of type 1 plaquettes that can be placed 
on a pyrochlore lattice is j. The 1 — 8 states do not realize 
this limit, but come fairly close to it with i of the plaque- 
ttes in the type 1 configuration. Therefore, the \f?> x \/3 
states are a "compromise" between the energy gain of 
Vi (which favors the type 1 plaquettes) and the energy 
loss of Vi (which disfavors the type 2 plaquettes), tak- 
ing into account the constraint x\ < X2- This becomes 
evident when calculating the energy of the various cell 
types. The state realizing the maximum fraction of type 
1 plaquettes found in RefjiS is comprised of only type 2 
cells. For s = |, type 1 cells are far more favorable in 
energy than type 2 cells, and type 1 cells are abundant 
in the 1 — 8 states. 
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